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ABSTRACT 

Aims. We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic BO - B5 supergiants of luminosity 
classes la, lb, lab, and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model 
predictions to observed UV spectral features is made. 

Methods. Fundamental parameters (e.g., T^g, logL,, mass-loss rates and CNO abundances) are derived for individual stars using 
CMFGEN, a nLTE, line-blanketed model atmosphere code. The impact of these newly derived parameters on the Galactic B supergiant 
Teff scale, mass discrepancy, and wind-momentum luminosity relation is examined. 

Results. The B supergiant temperature scale derived here shows a reduction of about 1 000 - 3 000 K compared to previous results 
using unblanketed codes. Mass-loss rate estimates are in good agreement with predicted theoretical values, and all of the 20 BO - 
B5 supergiants analysed show evidence of CNO processing. A mass discrepancy still exists between spectroscopic and evolutionary 
masses, with the largest discrepancy occurring at log (L/Lq) ~ 5.4. The observed WLR values calculated for BO - BO. 7 supergiants are 
higher than predicted values, whereas the reverse is true for Bl - B5 supergiants. This means that the discrepancy between observed 
and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass-loss rates as for O stars. The most surprising result 
is that, although CMFGEN succeeds in reproducing the optical stellar spectrum accurately, it fails to precisely reproduce key UV 
diagnostics, such as the N v and C iv P Cygni profiles. This problem arises because the models are not ionised enough and fail to 
reproduce the full extent of the observed absorption trough of the P Cygni profiles. 

Conclusions. Newly-derived fundamental parameters for early B supergiants are in good agreement with similar work in the field. 
The most significant discovery, however, is the failure of CMFGEN to predict the correct ionisation fraction for some ions. Such 
findings add further support to revising the current standard model of massive star winds, as our understanding of these winds is 
incomplete without a precise knowledge of the ionisation structure and distribution of clumping in the wind. 

Key words, stars: early type - stars: supergiants galaxies: Milky Way - stars: atmospheres stellar winds - stellar evolution - mass 
loss - abundances 



1. Introduction 

The study of luminous, massive stars is fundamental to im- 
proving our understanding of galactic evolution, since the 
radiatively driven winds of these stars have a tremendous 
impact on their host galaxies. This huge input of mechanical 
energy is responsible for creating H II regions, making a 
significant contribution to the integrated light of starburst 
galaxies and providing star formation diagnostics at both low 
and high redshifts. They substantially enrich the local ISM 
with the products of nucleosynthesis via their stellar winds and 
supernovae explosions and are a possible source of gamma ray 
bursts. It is therefore imperative to obtain accurate fundamental 
parameters for luminous massive stars since they contribute to 
many currently active areas of astrophysical research. 
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However there are still some uncertainties regarding the post- 
main sequence evolution of massive stars since their evolution 
is controlled by variable mass loss from the star as well as 
rotation, binarity and convective processes, the latter leading 
to surface enrichment as the products of nuclear burning are 
brought to the surface. Until recently, stellar evolution models 
failed to predict the correct amount of CNO processing in 
massive stars. However, new evolutionary tracks that account 
for the effects of rotation (" Maeder & Mevnel20 0r) show better 
agreement between predicted and observed amounts of CNO 
enrichment in massive stars. A far greater problem in stellar 
astrophysics is the determination of accurate observed mass loss 
rates. Recent research (e.g., Fullerton et al. 2006; Prinja et aU 
2005: M assa et al.ll2003l: IPuIs et a l. 2006; Repolust et alfeOOlT 
has shown that current OB star mass loss rates might be over- 
estimated by at least a factor of 10. Such large uncertainties 
in the mass loss rates of massive stars suggest that our under- 
standing of their winds is incomplete. It is now widely accepted 
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that the winds of both O and B stars are highly structured and 
clumped and therefore our assumptions that they are smooth 
and homogeneous are invalid. Evidence to support this claim 
has come from hydr odynamical, time-d e pendent simulations 
of ste llar winds (e.g. lOwocki et alj|l988t iRunacres & Owockil 
I2OO2I) . the latter of which proposed the idea that instabilites in 
the line-driving of the wind can produce small-scale, stochastic 
structure in the wind. Further evidence for the inhomogeneity 
of stellar wi nds comes in th e form of various o bservational 
studies (e.s.. iPuls et all 1200 6': 'Bour et et al.1 l2005t iMassa et all 
120031: iPrinia et al.ll2002l; iBianchi 200^ Time-series analyses of 
both Balmer and metal spectral lines (e.g., iPrinia et aLll2004l) 
in OB stars highlight clear, periodic patterns of vari abihty that 
corr espond to the evolution of structure in the wind. iPuls et aU 
( I2OO 6) demonstrated that the discrepancy between values of 
Mua and Mradio implies the presence of different amounts of 
clumping at the base of and further out in the wind. Massa et"aLl 
(12003,) showed that for a sample of O stars in the LMC, the 
empirical ionisation fractions derived were several orders of 
magnitude lower than expected, indicating a lack of dominant 
ions in the wind. Similar results were found by Prinja et al. 
(2005) for early B supergiants. More recently, iFuUerton et al.l 
(I2OO6) demonstrated that the ionisation fraction of Pv, which 
is dominant over a given range in T^ff in the O star spectral 
range, never approaches a value of unity. They subsequently 
showed that a reduction in mass loss rate of at least a factor of 
100 is required to resolve the situation. We intend to re-address 
the issue of the i onisation structure of early B supergiants, 
following on from IPrinia et alj (I2005h . in a forthcoming paper 
(Searle et al., 2007b, in preparation; hereafter Paper II). Such 
drastic reductions in OB star mass loss rates would have severe 
consequences for the post-main-sequence evolution of these 
stars; in particular it would affect the numbers of Wolf-Rayet 
stars produced and the ratio of neutron stars to black holes 
produced in the final stages of massive star evolution. 

Early type B supergiants are particularly important since they 
are the most numerous massive luminous stars and are ideal 
candidates for extra-galactic distance indicators, essential for 
calib rating the Wind -Momentum-Luminosity Relation (WLR) 
(e.g.. lKudritzki et aLl[T999.) . Research into this WLR calibration 
has highlighted a spectral type depende r ice for Galactic OBA 
type sta rs (e.g., [ Kudritzki et"an 1 19991: iRepolust etal] 12004 
[Markova et alj l20b4). whilst others have explored the effect of 
metallicity on the WLR by studying OB stars in the metal-poor 
enviro r iment of the Magell a nic Clouds (e .g.. iKudritzki & Pulsl 
I2OOOI; iTrundle et al.1 120041; lEvans et alj |2004b|). Accurately 
derived mass loss rates are essential in calibrating the WLR, 
yet discrepancies still exist between observed mass loss rates 
obtained from different wavelength regions (i.e. optical, UV 
or IR). Furthermore acknowledged discrepancies of up to 
30 % have been found bet ween observati onal and theoretically 
predicted mass loss rates ( Vink et a 012000). Vink has remarked 
that these discrepancies for the mass loss rates of B stars can 
be attributed to systematic errors in the methods employed to 
derive the observed values. Good agreement was fou nd between 
observed and predict ed mass loss rate s for O stars in IVink et alj 
(2000). Additionallv. iPuls et all (12006 ) recently derived values 
of both Miia and Mjadio, highlighting a discrepancy of roughly a 
factor of two between both values. 

The layout of this paper is as follows. ^ introduces the sam- 
ple of 20 Galactic BO - B5 supergiants upon which this analysis 
is based. |3] describes the methodology used in deriving funda- 



mental parameters for this sample. The results are presented in 
^ and a critical examination of the CMFGEN model fit to the 
UV spectra of the 20 B supergiants is made in ^ Finally the 
conclusions are given in ^ 

2. Observations 

Optical and UV spectra have been collected for a sample of 20 
Galactic B supergiants, covering the spectral range of B0-B5 
and including la, lb, lab and II luminosity classes as well as a 
hypergiant. Stars were only included in the sample if both opti- 
cal and lUE data were available for them. Where possible, stars 
were selected such that there would be 2 different luminosity 
classes at each spectral sub-type. The details of observational 
data for each star are given in Table 12.1] Fourteen of the chosen 
B supergiants belong to OB associations (Humphrevs 197^, 
so for these stars, the absolute magnitude given in Table 12.11 
is based on the distance to the relevant association; for the 
remaining six stars My and therefore the distance modulus is 
derived from photometry. 

For fifteen of the twenty B supergiants in our samp le, the optical 
spectra were taken from an existing data set (see ' Lennon et al.l 
1992 for further details). The blue spectra were observed using 
the 1-m Jacobus Kapteyn Telescope (JKT) at the Observatorio 
del Roque de los Muchacos, La Palma in October 1990 with the 
Richardson-Brealey Spectrograph and a R1200B grating. They 
have a wavelength coverage of 3950 - 4750 A, a spectral resolu- 
tion of 0.8 A and a signal-to-noise ratio ~ 150. The red spectra 
were obtained with the 2.5-m Isaac Newton Telescope (INT) us- 
ing the Intermediate Dispersion Spectrograph (IDS) and cover a 
wavelength range of 6260 - 6870 A ( iLennon etan fl992). The 
spectral resolution is 0.7 A and the signal-to-noise ratio is > 
100. Spectra for HD 192660, HD 185859, HD 190066 and HD 
191243 were taken at the INT in July 2003, again using the IDS. 
The R400B grating was used for the blue spectra, giving a cen- 
tral wavelength Ac of 4300 A, whereas the R600R grating was 
used for the red spectra, giving Ac ~ 6550 A. The signal-to-noise 
ratio was > 100 and the spectral resolution was 0.7 A. Finally 
high resolution time-average d blue and red spectra of HD 64760 
were provided by RKP (see lKaufer & Stahlll2002l for more de- 
tails). These spectra were taken in 1996 on the HEROS fiber- 
linked echelle spectrograph, which was mounted on the ESO 
50-cm telescope at the La Silla, Chile. The blue spectra have 
a range of 3450 - 5560 A whilst the red have 5820 - 8620 A. 
Both had a resolving power of 20000. The signal-to-noise ra- 
tio varied with lambda but for a red spectrum with a 40-minute 
exposure it was typically > 150. 

2.1. Sample details 

As previously mentioned, our sample covers a range of B0-B5 
supergiants with luminosity classes varying from la down to 
II in a cou ple of cases. Spectral type classifications have been 
taken from iLennon et al.l (Il992h . which includes some recent 
revisions. HD 204172 has been changed from BO lb to BO. 2 
la due to the strength of its Siiv lines and narrowness of its 
H lines. Comparing its UV resonance lines to those of the BO 
lb star HD 1 64402 supports t he change to a more luminous 
spectral type (Pri nia et alj|2002l) . Also both HD 164353 and HD 
191243 have been re classified from B5 la to B5 II stars. It is 
also worth noting that lde Zeeuw et af.l (1999) has questioned the 
membership of the stars HD 53138 and HD 58350 to ColHnder 
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Table 1. Observatio nal data for the sample of 20 Galactic B Supergiants. Spectral types and V magnitudes are taken from 
iLennon et al] ( 11992 ') for all stars except HD 192660, HD 18585 9, HP 190066 and HD 64760 The references for the spectral 
types of these 4 stars are as follows: HD 1 92660 from Walborn ( 197 lb : HD185859 from|Lesh| (Il968h : HD 190066 from Hiltne^ 
(119561) and HD 64760 from lHoffleit & Jas chek ( 1982) (from which the V magnitude of H D 64760 w as also taken). V magnitudes 
for the remaining 3 stars were obtained from Fernie ( 1983). Values of (B - V)o taken fro m Fitzgerald! ( 1970). Abso lute visual mag- 
nitudes, distance m oduli and cluster associations have been taken from: 1 . iBrown et~aL I (119941) . 21 iGarmanv & Sten cel ( 1992Dor3. 
iHumphrevi (Il978l) . For stars not associated with a cl uster, an absolute vis ual magnitude scale based on spectral type ( Egrea[T978l) 
was used. L1992 refers to archive data obtained from lLennon et alj (Il992l) . INT2003 denotes data taken on the 2.5m INT and RKP 
marks data supplied by R.K. Prinja. 



HD no. 


Alias 


Sp. Type 


V 


B-V 


Mv 


Association 


m-M 


Optical 


lUE 


37128 


e On 


BO la 


1.70 


-0.19 


-6.95 


Ori OB lb 


7.8' 


LI 992 


SWP30272 


192660 




BO lb 


7.38 


0.67 


-7.0 


Cyg 0B8 


11.8^ 


INT2003 


SWP44625 


204172 


69 Cvs 


B0.2 la 


5.94 


-0.08 


-6.2 


Cvs 0B4 


6.2 


LI 992 


SWP48900 


38771 


K Ori 


B0.5 la 


2.04 


-0.18 


-6.51 


OriOBlc 


8.0' 


L1992 


SWP30267 


185859 




B0.5 la 


6.48 


+0.44 


-6.6 




13.08 


INT2003 


SWP47509 


213087 


26 Cep 


B0.5 lb 


5.46 


+0.37 


-6.2 


Cep OBI 


12.7 


LI 992 


SWP02735 


64760 




B0.5 lb 


4.24 


-0.15 


-6.2 




10.44 


RKP 


SWP53781 


2905 


K Cas 


BC0.7Ia 


4.16 


+0.14 


-7.09 


Cas OB 14 


10.2^^ 


LI 992 


SWP54038 


13854 


V551 Per 


Bl lab(e) 


6.47 


+0.28 


-6.73 


Per OBI 


11.82 


LI 992 


SWP02737 


190066 




Bl lab(e) 


6.53 


+0.18 


-6.1 




12.63 


INT2003 


SWP18310 


190603 


V1768 Cyg 


B1.5 Ia+ 


5.64 


+0.54 


-7.0 




12.56 


LI 992 


SWP4325 


193183 




B1.5 lb 


7.01 


+0.44 


-6.24 


Cyg OBI 


11.3' 


LI 992 


SWP52716 


14818 


V554 Per 


B2 la 


6.25 


+0.30 


-6.93 


Per OBI 


11.8' 


L1992 


SWP 18658 


206165 


V337 Cep 


B2 lb 


4.74 


+.30 


-6.44 


Cep 0B2 


13.2 


LI 992 


SWP06336 


198478 


55 Cyg 


B2.5 la 


4.84 


+0.40 


-6.43 


Cyg 0B7 


9.6^ 


LI 992 


SWP38688 


42087 


3 Gem 


B2.5 lb 


5.75 


+0.22 


-6.26 


Gem OBI 


10.9' 


LI 992 


SWP08645 


53138 


24 CMa 


B3 la 


3.01 


-0.11 


-7.1 




10.11 


LI 992 


SWP30271 


58350 


T/CMa 


B5 la 


2.41 


-0.07 


-7.0 




9.41 


LI 992 


SWP30198 


164353 


67 0ph 


B5 Il(Ib) 


3.97 


+0.02 


-4.2 


Coll 359 


6.5 


LI 992 


SWP08560 


191243 




B5 Il(Ib) 


6.09 


+0.16 


-6.5 


Cyg OB3 


11.8'' 


INT2003 


SWP07737 



121 on account of insignificant proper motion and small 
parallax respectively. Seven of the twenty stars in our sampl e 
have been examined for Ha variability by iMorel et aP (l2004l) . 
who obtained both photometric and spectroscopic data on these 
objects in order to ascertain the am ount of variability p resent in 
their light-curves and Ha profiles. iMorel et akl (|2004|) quantify 
the amount of spectral and photometric variability observed as 
well as determining periods where cyclic behaviour is observed. 
The sample includes a rapid rotator, HD 64760 (discussed 
in detail in the following section) and a hypergiant, the B1.5 
Ia+ star HD 190603. Furthermore, there are several objects of 
interest in the sample for which a significant amount of re- 
search has akeady been undertaken and merit further discussion. 

- e Ori: There are several intriguing aspects of this star 
that are worth mentioning. Firstly i t has bee n noted as 
moderately nitrogen deficient by WalbornI (11976 ). Secondly 
it is known to un dergo significant variations in Ha, with 
iMorel et al] (|2004|) reporting variations of 81.9 %. Further 
studies by Prinja et al. (2004) have revealed variability in 
not only Ha but also H/3, He absorption and metal lines 
with a 1.9 day period. A modulating S-wave pattern has 
been discerned in the weaker lines, which cann ot be fully 
expla ined by current non-radial pulsation models (iTownsendl 
Il997h . These results highlight a direct connection between 
photospheric activity and perturbations in the stellar wind. 
Finally this star is the only normal early B supergiant to 
have a measured thermal radio flux (Blomme et al. 2002), 
from which a radio mass loss rate of log M = -5.72 was 



derived. 

- K Ori: Like e Ori, k Ori also exhibits spectra l variabihty in 
Ha. I ts Ha profile was studied in detail by iRusconi et al.l 
11980) who described it as a double-peaked absorption 
profile with a central emission core and broad profile 
wings, with variations on long (of the order of days) and 
short (of th e order of minutes) time scales. More recently 
iMorel etal] (|2004) reported chan ges in the p rofile amplitude 
and morphology of 32.6 %. WalbornI (11976 ) noted it as an 
example of a morphologically normal B supergiant in terms 
of the relative strengths of its CNO spectral lines. 

- HD 192660: 'Bidelman| (1988) noted that this star showed a 
"faint Ha emission with a slight P Cygni absorption", based 
on observations taken at the Lick observatory in 1957. Our 
spectrum of this star shows a very similar Ha profile that is 
also only weakly in emission. Walborn (1976) and Schill 
( 11985) both noted that HD 192660 displays evidence for 
nitrogen deficiency. 

- HD 64760: This star is classified as a rapid rotator, having 
a V sin / of 265 km/s. Many interesting studies have been 
carried out regarding the periodic an d sinusoidal mod- 
ulations of its optic a l and UV lines ('M assa et"aLl 119951: 
iFuUerton et all 119971: iKaufer & S tahl 2003), which have 
turned HD 64760 into a key object for improving our 
understanding of the spatia l structure an d varia tions of hot 
star winds. More recently, iKaufer et alj (l2006l) found, for 
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the first time, direct observational evidence for a connection 
between multi-periodic non-radial pulsations (NRPs) in the 
photosphere and spatially structured winds. More specifi- 
cally, they can use the interference of multiple photospheric 
pulsation modes on hourly timescales with wind modulation 
periods on time scales of several days. A beat period of 
6.8 days seen in the photosphere and base of the wind 
does not match with the derived periods of 1.2 and 2.4 
days for wind variability, being closer to the longer 5-11 
day repetitive timescales observed for discrete absorption 
components (DACs) in the lUE data sets. Evidently the 
precise nature of the wind-photosphere connection in this 
star is a complex o n e. Usi ng hydrodynamical simulations, 
ICranmer & Owockil (Il996h succeeded in confirming the 
existence of co-rotating interaction regions (CIRs). These 
are spiral structures in the wind that are produced through 
the collisions of fast and slow streams rooted in the stellar 
surface. 



K Cas: 'WalbornI ( 1 19721) classified this star as carbon rich, 
giving it a spectral type of BC0.7 la, when comparing its 
optical spectrum with that of HD 216411, a BO. 7 la star. He 
found that the nitrogen lines in k Cas were barely detectable, 
whereas the O II - C III blends were very prominent. 
Since Walborn also recognised that HD 216411 possesses 
a well developed nitrogen spectrum, he described k Cas 
as carbon-rich, rather than nitrogen-weak, with respect 
to a morphologically-normal B supergiant. Walborn had 
previously suggested that all OB stars begin their post- 
main-sequence evolution with an enhancement of carbon 
CWalborn 1971), which then becomes depleted as the star 
evolves and produces nitrogen as a by-product of the CNO 
bi-cycle. Therefore the implication of classifying k Cas as a 
BC 0.7 la star is that it is less evolved than other stars in the 
sample. Please consult Section 14.51 for a discussion of the 
carbon rich status of this star based on the results presented 
in that section. 



HD 13854: iMcErLean et al.1 (1 19991) describe this star as 
'highly processed' i.e. displaying a large amount of CNO 
processing in its spectrum. Its Ha profi le is seen m o stly in 
emission, assuming a P Cygni shape. iMorel et al.1 (12004 
found that not only does the Ha profile of this star vary 
by 47.3 %, but that these variations have a period of 1.047 
days +0.01. Hipparcos light curves also show evidence for 
periodic behaviour. 

HD 14818: possesses an Ho- profile with a P Cygni profile. 
iMorel et al. (2004) report variations of 34.8 % in the Ho- 
profile and, like HD 13854, they find that this star also 
shows periodic behaviour in its Hipparcos light curve. 

HD 42087: also has its Ha profile in emission but more 
importantly Mor el et al.l ((2004) reported significant Ha 
variability of 91.2 % {greater than the percentage variability 
that they found for e Ori), for which they find strong 
evidence of cyclic behaviour on a periodicity of 25 days 
+ 1-4 days. Moreover, the Ha variability correlates with 
variability in the He i 6678 A line, such that as Ha emission 
increases. He i 6678 becomes weaker. iMorel etaP (l2004l) 
also find periodic variability in its Hipparcos light curve. 



4x10" 



IxlO" 



Combined TLUSTY-CMFGEN density structure 
for HD 1B5B59 BO. 5 la 

Both models have the parameters 
Tetf = 36000K & log g = 3.13 
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log T„, 



Fig. 1. Temperature structure against Rosseland optical mean 
depth for the hydrostatic density structure produced from com- 
bining the subsonic TLUSTY velocity structure with the super- 
sonic CMFGEN velocity structure, such that the velocity and ve- 
locity gradient are constant. The above example is for the BO. 5 
la star HD 185859 {T^s^ 26 000 K, log g = 3.13) 

authors associating this star wi th OBN/OBC gr oups. It 
undergoes 66.4 % Ha variabihty dMorel et al.ll2004l) . 



3. Derivation of fundamental parameters 



Fundamental parameters were 
pie of st ars using the nLTE 



derived for this sam- 
stellar atmosphere codes 



HD 53138: IWalbornI (Il97(^ notes that this star shows 
a morphologically normal CNO spectrum, despite other 



TLUSTY dHu benv & Lmij [T99l lLan3 l2003b and CMFGEN 
(iHillier & Mil ler 1998). The application of TLUSTY, a plane- 
parallel photospheric code that does not account for the presence 
of a stellar wind, to modelling supergiants is valid as long as it is 
used solely to model purely photospheric lines. A n existing grid 
of B star TLUSTY models (iDufton et al.ll2005h was used as a 
base for this work and the grid (originally incremented in steps 
of ~ 2 000 K in Teff and 0.13 in log g) was refined further by 
RR in the range 15 000 K < T^s < 23 000 K to cover the same 
parameter space as the CMFGEN grid. The TLUSTY models 
provide a hydrostatic structure that can be input into CMFGEN, 
since the latter code does not solve for the momentum equation 
and therefore requires a density/velocity structure (see e.g., 
IHillier et al.ll200l iBouret etani200l iMartins et al.ll2005h . The 
TLUSTY input provides the subsonic velocity structure and 
the supersonic velocity structure in the CMFGEN model is 
described by a y6 - type law. The two structures are joined to a 
hydrostatic density structure at depth, such that the velocity and 
velocity gradient are consistent. The resulting structure is shown 
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in Fig. [T] which shows the change in Rosseland mean opacity, 
TRosj,with temperature and ensures that the model is calculating 
deep enough into the photosphere to sample the regions where 
the appropriate photospheric lines form (around -2 < logT«ojs 
< 0). This structure is then input into the CMFGEN model, 
adopting a yS-type velocity law of the form: 

vo + (vc.-vo)(l -R.lrf 

= — (1) 

1 + -!:2_e — 

where vq is the photospheric velocity, Vcore is the core velocity, 
Voo is the terminal velocity, /Zgff is the scale height expressed in 
terms of /?» and p is the acceleration parameter. The value of y6 
is normally determined from fitting the Ho' profile, as discussed 
in 93.21 assuming typical values of 1 .0 < yS < 1 .5 for B super- 
giant. Additionally, values of Vcore - 0.002 km/s and vq =0.1 
km/s are adopted for B supergiants as suggested by D.J. Hillier 
(priv. comm.). All parameters except log g were derived using 
CMFGEN and the precise details of the methods employed will 
be discussed next. The general method employed was to produce 
a grid of CMFGEN models of varying temperatures and lumi- 
nosities (all other parameters were kept constant), incremented 
in steps of 1 000 K in T^^ and 5 in log (L/Lq), and compare these 
synthetic spectra to observed spectra in order to constrain these 
two parameters ( 93.11 ). Once satisfactory values had been de- 
rived for the temperature and luminosity of a star, the mass loss 
rate, /3 velocity law, turbulent velocity (vtuib) ( 93.21 ) and CNO 
abundances ( 94.5b could be constrained. These parameters are 
all sensitive to changes in temperature and luminosity so can 
only be derived once those values have been fixed. Consistency 
checks are made after each parameter is altered to ensure that 
the resulting model output has not worsened the fit to the diag- 
nostic lines and overall spectrum. Given recent improvements in 
the treatment of metal lines in nLTE stellar atmosphere codes 
such as CMFGEN or FASTWIND, it is now possible to obtain 
accurate abundances by fitting (by eye) appropriate diagnostic 
lines for each element. All models were calculated to include 
the following elements; H, He, C, N, O, Al, Mg, Si, S, Ca & Fe, 
assuming solar abundances for silicon, magnesium, aluminium, 
phosphorous, sulphur, calcium and iron and adopting a relativ e 
number fraction of 5:1 for H:He (as used by Hillier et al. 2003|). 
CMFGEN adopts the superlevel approximation (see AndersonI 
119891 for more details), where a superlevel can consist of several 
or even many energy levels grouped together, such that all real 
levels j that form superlevel J have the same nLTE departure co- 
efficient (i.e., each component j is in Boltzmann equilibrium with 
respect to the other components). Details of the model atoms, 
including their full level and superlevel groupings, are given in 
Table 12 

3.1. Derivation ofT^s, log (L/Lq), iog g and CNO 
abundances 

In B stars, the silicon lines are used as the primary tempera- 
ture diagnostics, having the advantage that the abundance is well 
known as silicon is unaffected by nuclear processing. For BO - 
B2 supergiants the Si iv 4089 A and Si iii 4552, 4568 and 4575 A 
lines provide the main temperature diagnostics. Si iv 4089 A de- 
creases in strength as the temperature decreases until it is barely 
detectable at a spectral type of B2.5, which corresponds to T^ff 
~ 18 000 K. At this point the Si ii 4128-30 A doublet is present 
and replaces Si iv 4089 A as the main temperature diagnostic for 
B2.5 - B9 stars, along with Si iii 4552, 4568 and 4575 A. The 



Table 2. CMFGEN model atomic data, showing the number 
of full levels and superlevels treated as well as the number of 
bound-bound transitions considered for each ion included in a 
CMFGEN model. 



Ion 


Full Levels 


Superlevels 


b-b transit 


HI 


30 


20 


435 


He I 


59 


41 


590 


Hell 


30 


20 


435 


C II 


53 


30 


323 


C III 


54 


29 


268 


C IV* 


18 


13 


76 


NI 


22 


10 


59 


Nil 


41 


21 


144 


NIII* 


70 


34 


430 


O I 


75 


18 


450 


O II 


63 


22 


444 


III* 


45 


25 


182 


M<y II 


45 


18 


362 


Al II 


44 


26 


171 


Al III 


65 


21 


1452 


Sill 


62 


23 


365 


Sim 


45 


25 


172 


SilV 


12 


8 


26 


SII 


87 


27 


786 


SIII 


41 


21 


177 


SIV* 


92 


37 


708 


Call 


12 


7 


28 


Fell 


510 


100 


7501 


Felll 


607 


65 


5482 


FelV 


272 


48 


3113 


Fe V* 


182 


46 


1781 



He I lines at 4144 A, 4387 A, 4471 A and 4713 A and Mgii 
line at 448 1 A can also be used as secondary criteria for both 
temperature and luminosity, since they are sensitive to changes 
in both parameters. The principal luminosity criteria used in 
spectral classification of B stars is the ratio of Si iv 4089 A to 
He I 4026 A, 4121 A and/or 4144 Afor BO - Bl supergiants, 
whereas for stars later than Bl the ratio of Si iii 4552, 4568 and 
4575 A to He I 4387 A is used. The procedure adopted for de- 
riving values of T^g, log (L/Lq), log g and CN O abundances is 
the same method adopted by Hillier et al. (2003) jCrowther et al.l 
(l2006h : lBouret et all (l2003h : IMartins et alj (i20(j5h and is as fol- 
lows: 

1 . An optical stellar spectrum of a chosen star is compared to a 
grid of CMFGEN models which differ in values of Teff and 
luminosity (all other parameters are kept constant). 

2. A value of Teff is selected for the star by finding the model 
that provides the best fit to the temperature sensitive silicon 
lines. The diagnostic lines are Siiv 4089 and Sim 4552, 
4568 and 4575 for B0-B2 supergiants and Sin 4128-30 and 
Sim 4552, 4568 and 4575 for B2.5-B5 supergiants. This 
temperature would then be confirmed by checking that the 
model also succeeded in fitting the helium lines listed earlier 
and Mg ii (but these lines were not used to derive the initial 
value of Teff) to ensure it provided a reasonable match to all 
temperature-sensitive lines. 

3. Once a value of Teff has been chosen, an inital estimate of 
log (L/Lq) is made by selecting the model from the grid 
(at the chosen value of Teff) whose value of My (output 
from the model) best matches the observed value of My 
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for the star in question. The luminosity is then constrained 
further by taking observed values of My, the absolute visual 
magnitude and V, together with the estimate of A(V), were 
used to obtain an initial estimate of the distance modulus. 
Optical photometry and ultraviolet spectroscopy were then 
de-reddened with respect to the model spectral energy 
distribution to obtain revised estimates of E(B-V) and My- 
This My derived from the model was then compared to an 
observed My, if the values matched then the model lumi- 
nosity was correct. If not, the value of My was translated 
into a bolometric correction to obtain an estimate of the 
corrected luminosity for the model, which was then rerun 
with this value for the luminosity. This iterative process was 
continued until the observed and model values of My were 
in reasonable agreement and then the resulting model was 
checked against the T^ff diagnostic lines (Sin, Sim, Siiv) 
and derived value of T^ff was adjusted if necessary. 

4. Estimates of log g were then made from TLUSTY fits to the 
Hy lines of the observed spectra. Hy is normally chosen as 
the log g diagnostic since Ho' and H/3 suffer from too much 
wind fill emission; H6 was used as a secondary diagnostic 
to check for consistency with values of log g derived from 
Hy. Since Hy is affected by an O ii blend around 4350 A and 
H(5 has a blend with N iii 4097, neither line is ideal but both 
implied the same log g values so we can have confidence 
that the derived log g values are not affected by these blends. 
The adopted log g value was then incorporated into the 
CMFGEN model and again the derived T^s and log (L/Lq) 
values were revised if necessary. 

5. Next CNO abundances were derived by varying the abun- 
dance of each element until the appropriate diagnostic lines 
were fitted by the model. For nitrogen, we use the N iii 4097 
A line for BO - B2 stars (which is blended with H6) and the 
Nil 3995 A, Nil 4447 A and 4630 A lines (for all B stars) 
as primary diagnostics. The main diagnostics for carbon and 
oxygen are the 4267 A and 4367 A lines respectively; how- 
ever their abundances are confirmed by checking for good 
fits to On blends at 4070 A, 4317 - 4319 A and 4650 A, the 
On fines at 4590 A, 4596 A and 4661 A and the Cii dou- 
blet at 6578, 6582 A. The errors on constraining the CNO 
abundances using this method were typically up to ~ 0.3 dex. 



1. Values for the mass loss rate of each star were constrained 
by fits to the Ha profile, with each fit aiming to reproduce 
the overall profile shape and amplitude. 

2. Values of Voa are best determined from the UV s o values ob- 
tained through UV line synthesis modelling (see lPrinia et al.l 
I2OOI Paper II) are used here. 

3. The value of /3 is also varied in order to improve the shape 
of the model profile fit with respect to the observed profile, 
but this has no effect for Ha profiles in absorption, in which 
case values obtained from SEI analysis were used. 

4. Estimates of the microturbulent velocity, Vmrb, were then 
made by fitting the Sim 4552, 4568 and 4575 fines. 

5. The model fit to the temperature and luminosity sensitive 
lines was checked after altering these parameters to ensure 
consistency. 



The sample of 20 B supergiants showed some variation in 
overall profile shape, with four stars showing Ha in emission, 
six in absorption and the remainder displaying a more complex 
morphology. In the last case, the profiles are partly in absorption 
with some emission component also detectable, implying that 
the profile has been filled in by stellar wind emission. The 
B2 - B5 supergiants display Ha profiles with a P Cygni profile 
shape, which implies that line scattering is playing a significant 
role in the line's formation, though this is not normally observed 
until late B/A supergiants. The majority of stars in this sample 
have no record of Ha variability so that any changes in the line 
profile morphology and amplitude can be considered negligible 
for the purpose of this analysis. However e Ori, HD 13854, 
HD 14818 and H D 42087 afi display significant Ha variabifity 
according to Mor el et al.| ^2004). In view of this problem, 
for £ Ori we have assumed a radio mass loss rate, M r adio, o f 
1.9 X lO"*" M0 yr"', as measured by iBlomme et all ( |2002|) . 
thereby avoiding the inaccuracies involved in deriving M from 
a variable Ha line. An estimate of the error in fitting the Ha 
profile of this star is given nonetheless is Table [8] Unfortunately 
this approach is not possible for the other 3 stars since there are 
no reliable M,adio values available in the literature. The problems 
associated with deriving M from the Ha profiles of these stars 
are discussed in the next section. 



3.2. Determination of Stellar Wind Properties 

The stellar wind parameters M, p and Vturbwere then con- 
strained using the usual method outlined below (again the 
same procedure use d by Hillier et al. 2003; Crowther et al.l2066t 
iBouret et al.l 120031: iMartins et al.i 120051) . CMFGEN allows for 
a treatment of turbulence in the stellar wind by assuming a 
radially-dependent microturbulent velocity, defined as 

(Vmax - Vniin)v'('-) 
Vtuib = Vmin + (2) 

Voo 

where Vn,;,, is the minimum turbulent velocity occurring in 
the photosphere a nd Vmax is the maximum turbulent velocity. 
iHilfieretal.! (120031) found that varying the turbulent velocity had 
little effect on the temperature structure calculated by CMFGEN. 
In this work, values of Vmin = 10 km/s and v^ax = 50 km/s are 
adopted as limits. Typically v reaches the value of Vturb around 
r = 1.057?.. 



4. Results 

4. 1 . The B supergiant Teff scale 

Model fits to the optical spectra of HD 192660 (BO lb), HD 
213087 (B0.5 la), HD 13854 (Bl lab), HD 193183 (B1.5 lb), 
HD 14818 (B2 la), HD 198478 (B2.5 la), HD 53138 (B3 la) 
and HD 58350 (B5 la) are shown in Fig. |2] (4050 - 4250 A), 
Fig. [3] (4250 - 4450 A) and Fig. g] (4450 - 4650 A). Overall 
CMFGEN has succeeded in providing excellent fits to the 
observed spectrum of each star The models succeed in repro- 
ducing the H, He, Si and Mg lines quite accurately. However, 
some individual spectral lines are more difficult to model than 
others. The Siiv 4089 A line is sometimes underestimated in 
BO - B2 supergiants, with the effect being most pronounced 
in Bl - B2 supergiants, which might be partly due to a blend 
with O II. It is also noticeable that the model Si iv line displays 
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Fig. 2. Overall CMFGEN fit to the optical spectrum of BO - B5 supergiants (4050 A - 4250 A). The solid black line is the observed 
spectrum and the red line denotes the CMFGEN model fit. 



a slight sensitivity to mass loss. iHilher et al.l (l2003l) also noted 
that some model photosphe ric lines c an be affected by mass loss 
and more importantly, Dufton et al.l (QQ05) noted when using 
FASTWIND that Siiv 4116 A and the Sim multiplet 4552, 
4568, 4575 A were affected by the stellar wind. However, whilst 
using CMFGEN, we have not observed any significant stellar 
wind effects on the Si ni multiplet. The values of T^ff derived for 
these stars can still be justified since the model spectrum still fits 
the rest of the spectrum, including the Si iii, Mg ii and He i lines, 
very well. In the cases where the model does underestimate the 
Si IV 4089 A line, the use of a model with a higher T^ff that 
provided a better fit to the Si iv line would provide a worse fit 
to the rest of the observed spectrum. Note that the values of 
Teff obtained in these cases were still derived using the silicon 
ionisation balance and the effect of the compromise attained 
between fitting the Si iv line underestimated by the model and 
the rest of the spectrum is reflected in the value of AFeff quoted 
in Table |3] It is also intriguing to note that CMFGEN predicts 
two absorption lines at 4163 A (see in B1.5 - B5 supergiants) 
and 4168.5 A (seen in all BO - B5 supergiants) that are not 



observed in any of the sample stars. These predicted hne s 
also appear in the CMFGEN models of ICrowther et all (|2006|) . 
where it appears that they have identified the line at 4168.5 A as 
He I but no explanation is given for the line at 4163 A. We can 
confirm the identify of the line at 4168.5 A and also add that the 
line at 4163 A is Fe m. 



The values of Tgff, log g, log (L/Lq), E(B-V) and My derived 
for each of the 20 B supergiants in the sample are listed in Table 
|3] These results show that BO - B5 supergiants have a range 
in TefF of 14500 - 30000 K, in log (L/Lq) of 4.30 - 5.74 and 
that their stellar radii vary from about 20-71 Rq. They also 
exhibit a range of -6.04 < My < -7.26 in brightness, con- 
firming their status as some of the brightest stars in our Galaxy. 
The temperature scale for B supergiants derived here is shown 
in Fig.|5] plotted against spectral type. A drop of up to 10000 K 
in temperature is witnessed between BO - Bl, whereas at lower 
spectral types, the Tetf scale shows a more grad ual decrease in 
reff- T he Galactic O star Teg scale published by iRepolust et al.l 
(12004 ranges from an 02 If star with T^g = 42 500 K down to 
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Fig. 3. Overall CMFGEN fit to the optical spectrum of BO - B5 supergiants (4200 A - 4450 A). The solid black line is the observed 
spectrum and the red line denotes the CMFGEN model fit. 



an 09.5 la star with = 29000 K and an 09.5 lb star with 
T^s- 32000 K, meaning that the B supergiant Teff scale pre- 
sented here carries on smoothly from the Galactic O supergiant 
Jeff scale. Similarly the B supergiant Tetj scale ends with B5 
lb stars having T^-s ^ 15 000 K and the Galactic A supergiant 
Teff scale derived by I Vend (119951) begins with a T^ff of 9950 K. 
A gap between the B and A supergiant Tgff scales is expected 
since none of the recently published B star T^ff scales include 
B6-9 stars. The B supergiant temperature scale derived here also 
demonstrates the difference in Tetf between B la and B Ib/Iab/II 
stars. BO - B2 Ib/Iab stars are found to be up to 2500 K hot- 
ter than BO - B2 la stars, with the exception of the stars HD 
190603 (B1.5 Ia+) and HD 193183 (B1.5 lb). However, a less 
significant difference of 500 K in T^ff is found between B2 - 
B5 la and B2 - B5 Ib/II stars, with the B lb stars again being 
hotter than their more luminous counterparts; this discrepancy 
is well within the margins of eiTor in deriving T^ff as typically 
ArefF= 500 - 1 000 K. We have compared our Galactic B super- 
giant Tgg scale to other published values (Trundle et al. 20(3 
Trundle & Lennon 2005; McErLean et al. 1999; CrowtheretaD 
2006 ) in Table |4] Where each author has several stars with the 
same spectral type, the values are averaged and marked with an 



asteris k in the table. Note that the results of iMcErLean et all 
(Il999h were obtained with an unblanketed stellar atmosphere 
code. If w e compare our derived T eft values with those of the un- 
blanketed |McE^eanetal]([i999|) Tsff scale, the use of a stellar- 
atmosphere code with a full treatment of line blanketing has 
the effect of lowering T^ff by 1 000 - 3 000 K for Galactic B 
supergiants. This is not as drastic as the reduction found for 
O super giants, which can be as high as 7 000 K for extreme 
stars (ICr owther et al.ll2002h. If we compare our derived Teff 's to 
those of McE rLean et al.l (1 19991) . with whom we have 10 target 
stars in common (HD 37128, HD 38771, HD 2905, HD 13854, 
HD 193183, HD 14818, HD 206165, HD 42087 and HD 53138), 
we find r easonably good ag reement except for B 1 la/Iabs, where 
the McE rLean et al.l d 19991) results imply that a Bl supergiant 
is 2500 - 3 000 K hotte r than our va lues. The SMC B super- 
giant temperature scale (iTrundle et al. 2004; T rundle & LennorJ 
2005) also implies a much hotter Bl supergiant, but it is ex- 
pected that SMC stars will be h otter than G alactic stars (see e.g. 
the O star temperature scales of lMassev et a l. (2005) (SMC) and 
'Repolust etaLl (|2004|) (Galactic) where the SMC stars are up to 
4 000 K hotter than the Galactic ones). 
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Fig. 4. Overall CMFGEN fit to the optical spectra of BO - B5 supergiants (4450 A - 4650 A). The soUd black line is the observed 
spectrum and the red line denotes the CMFGEN model fit. 



4.2. Log g estimates 

An example of a TLUSTY log g fit to the Hy profile of HD 
164353 is shown in Fig. [T] Fits to the Hy and Hd profiles 
of all 20 B supergiants can be found online. TLUSTY is a 
purely photospheric code therefore the model Balmer lines 
are also photospheric; however in B supergiants the observed 
Balmer lines suffer from wind contamination, where hydrogen 
photons emitted in the wind at the same wavelengths as Hy 
and H6 'fill in' the absorption profile. This makes it appear 
more 'shallow' when compared to a photospheric profile and 
explains the difference in depth between the two profiles shown 
in Fig. |2] Values of log g derived from Hy and H6 were in 
good agreement in general and no discrepancies larger than the 
margin of error on log g were found. Only small discrepancies 
were found for BO - BO. 5 stars where a model with a log g 
value 0.13 dex higher than the adopted value might provide 
a slight better fit to the H6 profile. However this ambiguity 
can be attributed to the influence of a large Niii blend on the 
blue wing of H6 masking where the actual wing of the profile 
should really lie. In these cases a very good fit is made to Hy 



so the value derived from Hy is taken. Some difficulties were 
encountered when trying to fit the Hy and H6 profiles of HD 
192660, HD 64760, HD 190603, HD 13854 & HD 190066 due 
to the observed asymmetry of the Hy and Hd profiles. This is 
particularly evident in the HD 190603, a B 1 .5 la hypergiant with 
a strong wind evident from the P Cygni shape of its Hfi profile. 
The Teff - log g scale derived from this work is shown in Fig. 
|6] where higher log g values are found for B lb stars. The log g 
values derived for B la stars are 0. 1 - 0.2 dex higher than t hose 
obtained by Kudritzki et al] (Il999h : ICrowther e t aT ('20 06|) for 
a samp l e of Galactic B supergiant s, whereas the Tr undle et al.l 
( 120041) : iTrundle & LennonI ( l2005h values for SMC B super- 
giants are generally higher than those for Galactic B supergiants. 



4.3. Stellar Masses and the Mass Discrepancy 

Using our estimates of log g, spectroscopic masses have been 
derived for each of the 20 B supergiants and imply a range 
of 8 < Mt < 52. Estimates of the evolutionary mass, Mevoi, 
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Table 3. Fundamental parameters (Teff, log ^, log (L/Lq), RtiRo ), E(B-V) Myand Mevoi) derived for the sample of 20 Galactic B 
supergiants. Values of sin / are taken from lHowarth et al] d 19971) and are expressed in km/s. 



HD no. 


Sp. type 


Teff (K) 




log (L/Lo) 




E(B-V) 


Mv 


Ve sin / 


Mevol 


37128 


BO la 


27500 ± 1000 


3.13 


5.73 + 0.11 


32.4 ± 0.75 


0.08 ± 0.02 


-6.89 ± 0.05 


91 


40 


192660 


BO lb 


30000 + 1000 


3.25 


5.74 + 0.13 


23.4 ± 1.03 


0.80 + 0.10 


-6.66 + 0.10 


94 


33 


204172 


B0.2 la 


28500 + 1000 


3.13 


5.48 + 0.27 


22.4 + 3.23 


0.12 ± 0.04 


-6.07 ± 0.30 


87 


41 


38771 


B0.5Ia 


26000 + 1000 


3.00 


5.48 + 0.14 


27.0+ 1.24 


0.07 + 0.01 


-6.48+0.10 


91 


33 


185859 


B0.51a 


26000 + 1000 


3.13 


5.54 + 0.14 


29.1 + 1.34 


0.53 + 0.02 


-6.54 + 0.10 


74 


35 


213087 


B0.5Ib 


27000 + 1000 


3.13 


5.69 + 0.11 


32.0 + 0.01 


1.30 + 0.01 


-6.20 + 0.10 


88 


40 


64760 


B0.5Ib 


28000 + 2000 


3.38 


5.48 + 0.26 


23.3 + 2.15 


0.15 + 0.05 


-6.36 + 0.20 


265 


33 


2905 


BC0.7Ia 


23500 + 1500 


2.75 


5.48 + 0.22 


33.0+ 1.52 


0.29 + 0.03 


-7.00 + 0.10 


91 


33 


13854 


Bl lab 


20000 + 2000 


2.50 


5.54 + 0.57 


49.2 + 1.33 


0.60 + 0.10 


-6.41 + 0.50 


97 


33 


190066 


Bl lab 


21000+ 1000 


2.88 


5.54 + 0.20 


41.4+ 1.89 


0.55 + 0.02 


-6.04 + 0.10 


82 


33 


190603 


B1.5 la+ 


19500 + 1000 


2.38 


5.41 + 0.23 


44.5 + 3.07 


0.70 + 0.05 


-6.85 + 0.15 


79 


27 


193183 


B1.51b 


18500+ 1000 


2.63 


5.00 + 0.26 


30.8 + 2.84 


0.70 + 0.06 


-6.43 + 0.20 


68 


18 


14818 


B2 la 


18000 + 500 


2.38 


5.40 + 0.27 


51.4 + 7.10 


0.62 + 0.10 


-6.70 + 0.30 


82 


26 


206165 


B21b 


18000 + 500 


2.50 


5.18 + 0.26 


39.8 + 5.50 


0.56 + 0.10 


-6.64 + 0.30 


73 


21 


198478 


B2.5 la 


17500 + 500 


2.25 


5.26 + 0.14 


46.1 +2.12 


0.40 + 0.05 


-7.26 + 0.10 


61 


23 


42087 


B2.5 lb 


18000+ 1000 


2.50 


5.11 +0.24 


36.6+1.69 


0.60 + 0.02 


-6.11 +0.10 


71 


21 


53138 


B31a 


16500 + 500 


2.25 


5.30 + 0.27 


54.7 + 7.56 


0.35 + 0.10 


-6.79 + 0.30 


58 


23 


58350 


B5 la 


15000 + 500 


2.13 


5.18 + 0.17 


57.3 + 2.64 


0.05 + 0.04 


-7.12 + 0.10 


50 


21 


164353 


B5Ib 


15500 + 1000 


2.75 


4.30+1.30 


19.6 + 8.05 


0.71+0.05 


-6.15 + 2.00 


44 


10 


191243 


B5Ib 


14500 + 1000 


2.75 


5.30 + 0.37 


70.8 + 3.26 


0.93 + 0.03 


-6.59 + 0.10 


38 


21 
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Fig. 5. The Galactic B supergiant T^g scale as a function of spec- 
tral type. B la stars are indicated by an asterisk, whilst B lb stars 
are marked by a plus sign. 



were then obtained using ou r derived stellar parameter s and the 
stellar evolutionary tracks of lMevnet & Maedeij(l2000b . The po- 
sitions of our 20 Galactic B supergiants on the Hertzsprung- 
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Fig. 6. The Galactic B supergiant T^ff - log g scale. 



Russell diagram, alo ng with other Gal actic B supergiants 
("Crowther et al.'2006), SMC B supergiants (Trundle et al. 2004t 
Trundle & Lennon 2005) and Galactic O stars (Repolust et al] 
2004), are shown in Fig. [8] Here, the iMevnet & 'Maeded (l2000h 
stellar evolutionary tracks have been used, which include the 
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Table 4. Values of T^ff (expressed in te rms of 10^ 
K) obtained in this thesis work and from Trundle et alj 
(12004 ); Trundle & Lennon (2005); Kudritzki et al. (199^ 
iMcErLean et al.l d 199911 . Values marked with an asterisk denote 
where values from one author have been averaged and are quoted 
to 1 decimal place, f the BO. 5 lb star HD 64760 has been omitted 
here because it is a rapid rotator. 



1.0 \ 



Sp type 


This work 


Trundle 


McErlean 


Crow 


BO la 


27.5 


27.0* 


28.5 


27.4* 


BO lb 


30.0 


- 


- 


- 


B0.2 la 


28.5 


- 


28.5 


- 


BO. 5 la 


26.0 


27.3* 


27.5 


26.0* 


B0.5 lb 


27.0t 




26.5* 




B0.7 la 


23.5 




24.0 


22.9 


Bl la 




23.8* 




22.0 


Bl lab/lb 


20.5* 




23.3* 


21.8* 


B1.5 la 


19.5 


21.3* 


21.25 


18.17 


B1.5 lb 


18.5 




22.3* 




B21a 


18.0 


19.0 


19.83 


18.6* 


B21b 


18.0 




20.8* 




B2.5 la 


17.5 


16.5 


18.0 


16.5 


B2.5 lb 


18.0 




20.5 




B3 la 


16.5 


14.0 


17.9* 


15.8* 


B4 lab 






16.5 




B5 la 


15.0 


14.5* 


15.4* 




B5 Ib/Il 


15.0* 




15.8* 





effects of rotation and are therefore more appropriate for OB 
supergiants. In order to demonstrate the effect of different stel- 
lar parameters on a star's precise position on the HR diagram. 
Galactic B supergiants common to both our sample and that of 
ICrowther et al. (2006) are joined by a dotted line. A compar- 
ison of both masses is shown in Fig. |9] For 14 out of the 2 
B supergiants, Mevoi > ^spec as found by iHerrero et alJ (l2002h . 
However, for the 5 other stars, which (excluding the rapid ro- 
tator HD 64760) have log (L/Lq) > 5.54, Mgvoi < M^pec- The 
dependence of the mass discrepancy with luminosity is exam- 
ined further in Fig. [TO] and compared to the mass disc r epanc y 
for SMC B s upergiants investigated by i Trundle et"an ( |2004|) ; 
iTrundle & Le nnon (2005). Both data sets exhibit a peak in the 
mass discrepancy at 5.4 < log (L/Lq) < 5.5 that drops off quite 
rapidly. 



4.4. Calibration of B supergiant fundamental parameters 

A calibration of stellar atmosphere parameters (i.e., Teff, logL,, 
log g, Mt, R,) according to spectral type has been carried out, 
using the fundamental parameters derived for our sample and 
that of lCrowther et al.l (12006) . A linear regression was applied to 
the trend of T^ff with spectral type; once the Teff scale had been 
established, linear regressions were made to the trends of log Teff 
vs. log Lt and log T^ff vs. log g, from which the values of and 
M, were then calculated. The resulting values of Teff, log L,, log 
g, Mt and for each spectral type are shown in Table |5] 



4.5. Evidence for CNO processing in B supergiants 

As previously mentioned, it was I WalbornI (Il976h who first 
suggested that the nitrogen and carbon anomalies found in OB 
stars can be explained by their evolutionary status, with OBC 
stars being the least evolved. It therefore follows that a typical 
OB supergiant should display some partial CNO processing, in 
the form of nitrogen enrichment accompanied by CO depletion. 



0.8 
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0.4 - 




0.2 I I I I I I I I I I I I I I I I I I I I I 
4334 4336 4338 4340 4342 4344 4346 
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Fig. 7. Example of a TLUSTY log g fit to Hy profile of the B5 
Ib/II star HD 164353. A value of log g = 2.75 is used here. 



Several authors ( iTrundle et alj |2004l; iTrundle & Lennonl l2005l; 
lEvans et alj|2004ar Vend 1995 ) have found N enrichments and 
CO depletions in OB A stars with respect to solar abundances. 
All 20 B supergiants in our sample show evidence for partial 
CNO processing in their spectra. The details of the CNO 
abundances derived for individual stars are given in Table|6] 

The majority of Galactic B supergiants show a modest nitro- 
gen enrichment, but some stars (HD 37128, HD 192660 and 
HD 191243) are slightly nitrogen deficient. Wal bornI (Il976h 
observed that Orion belt stars such as HD 37128 (e Ori) are 
nitrogen deficient due to the weakness of the N iii 4097 A and 
4640 A (blend) spectral lines. The largest nitrogen enhance- 
ments are seen for Bl - B2 stars (HD 13854, HP 190 603 and 
HD 14818). It is of interest to note that I WalbornI d 19761) classed 
HD 13854 as a morphologically normal B supergiant (as well as 
HD 38771), whereas we have found a modest yet significant N 
enrichment in this star 

In general, the fits to the CNO diagnostic Unes are good. On the 
whole, nitrogen abundances constrained from Nii 3995 A and 
N III 4097 A are in good agreement; the only exceptions being 
HD 37128 and HD 193183. It is well-documented in the litera- 
ture that there exists a discrepancy between carbon abundances 
derived from the Cii 4267.02, 4267.27 A multiplet line and the 
C II multiplets at 6578 and 6582 A, due to their strong sensitivity 
to nLTE effects and the adopted stellar parameters (see e.g., 
Nieva & Przvbilla 2006). A combination of high-resolution and 
high signal-to-noise spectra, along with sufficiently- detailed 
model atoms, are required to attempt to resolve this problem. 
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Table 6. Derived CNO abundances for the sample of Galactic B supergiants (expressed as log + 12). The amount of nitroge 
enrichment relative to carbon and oxygen respectively is given in the last two columns, calculated as log j ~ l^S (i^) ■ 



HD no. 


Sp. type 




Vf sin i 


C 


N 


O 


N/C 


N/0 


SUN 


G2 V 


5770 


- 


8.39 


7.78 


8.66 


-0.61 


-0.88 


37128 


BO la 


27500 


91 


7.66 


7.31 


8.68 


+0.26 


-0.49 


192660 


BO lb 


30000 


94 


8.02 


7.51 


8.73 


+0.10 


-0.40 


204172 


B0.2 la 


28500 


87 


7.66 


7.71 


8.66 


+0.66 


-0.07 


38771 


B0.5 la 


26000 


91 


7.74 


8.15 


8.73 


+ 1.02 


+0.30 


185859 


B0.5 la 


26000 


74 


7.72 


7.95 


8.53 


+0.84 


+0.30 


213087 


B0.5 lb 


27000 


88 


8.00 


8.15 


8.73 


+0.76 


+0.30 


64760 


B0.5 lb 


28000 


265 


7.99 


8.15 


8.73 


+0.77 


+0.30 


2905 


BC0.7 la 


23500 


91 


7.99 


8.16 


8.80 


+0.78 


+0.24 


13854 


Bl lab(e) 


20000 


97 


7.66 


8.51 


8.80 


+ 1.46 


+0.59 


1 9flO(S6 

1 "\J\J\J\J 


Rl Tahl'p'l 


21000 


82 


7.88 


8.15 


8.53 


+0.88 


+0.50 


190603 


B1.5 la+ 


19500 


79 


7.48 


8.76 


8.73 


+ 1.89 


+0.91 


193183 


B1.5 lb 


18500 


68 


7.66 


8.15 


8.73 


+ 1.10 


+0.30 


14818 


B21a 


18000 


82 


7.66 


8.72 


8.90 


+ 1.67 


+0.70 


206165 


B21b 


18000 


73 


7.96 


8.15 


8.43 


+0.80 


+0.60 


198478 


B2.5 la 


17500 


61 


7.86 


8.29 


8.45 


+ 1.04 


+0.72 


42087 


B2.5 lb 


18000 


71 


7.76 


8.11 


8.80 


+0.96 


+0.19 


53138 


B3 la 


16500 


58 


7.78 


8.32 


8.60 


+ 1.15 


+0.60 


58350 


B5 la 


15000 


50 


7.78 


8.29 


8.75 


+ 1.12 


+0.42 


164353 


B5 lb/11 


15500 


44 


7.78 


7.89 


8.53 


+0.72 


+0.24 


191243 


B5 lb/11 


14500 


38 


7.70 


7.65 




+0.56 





Table 7. Comparison of mean published CNO abundances for OBA supergiants (expressed as log j + 12). 



Author 


Stellar group 


C 


N 





N/C 


N/0 


This work 


Gal BSGs 


7.79 


8.12 


8.68 


+0.57 


+0.26 


Crowther et al. (2006) 


Gal BSGs 


7.87 


8.33 


8.47 


+0.68 


+0.59 


Trundle et al. (2004-5) 


SMC BSGs 


7.27 


7.71 


8.13 


+ 1.38 


1.04 


Evans et al. (2004) 


LMC OBSG 


7.49 


8.55 


8.02 


+ 1.28 


+ 1.26 


Evans et al. (2004) 


SMC OBSG 


7.30 


7.94 


8.01 


+0.86 


+0.66 


Venn(1995) 


Gal ASGs 


8.14 


8.05 




+0.13 




Venn(1999) 


SMC ASGs 




7.33 


8.14 




-0.08 



It is unlikely that our data is of a suitable resolution and 
signal-to-noise to attempt to solve this discrepancy, but we 
will nonetheless discuss our findings as appropriate. For our 
sample of B supergiants, the C ii multiplets at 6578 and 6582 
A are not prominent for BO - Bl supergiants; however for 
Bl - B5 stars the lines are distinguishable. The fits to the Cii 
4267.02, 4267.27 A multiplet are very good but CMFGEN 
tends to overestimate the Cn 6578 and 6582 A multiplets. In 
the case of constraining the nitrogen abundances, in general 
good agreement is found between the abundance implied from 
the Niii 4097 and Nii 3995, 4447, 4630 A. Some exceptions 
are found for some BO - BO. 5 supergiants (HD 37128, HD 
204172, HD 38771 HD 192660) where very good fits are made 
to Nil 3995 and 4447 A, but Niii 4097 and Nii 4630 are both 
underestimated, with the model producing a much weaker N ii 
4447 line than observed. Evidently increasing the nitrogen 
abundance would then cause the N ii 3995 and 4447 Alines to 
be overestimated. For HD 185859 and HD 213087, much better 
agreement is found between all four nitrogen diagnostics. 

There is still an intriguing contradict ion that ^ Cas, which has 
been defined as a carbon-rich star dWalbornl 1 1976 1) has very 
similar CNO abundances to the stars HD 64760, HD 213087 
which have not been noted as carbon rich by any other authors. 
The original criteria for classifying k Cas as a carbon-rich star 



were based on the weakness of its nitrogen lines as w ell as the 
strength of its carbon lines; this makes sense since (as IWalbornI 
1976 explains) it is expected that nitrogen deficiency will be 
accompanied by carbon enrichment. In order to resolve this 
discrepancy, the lUE spectrum of k Cas has been compared 
to the lUE spectrum of the B0.7 la star HD 154090 (see Fig. 
fTTT i. Looking at the Cn 1324 A line, it is certainly no stronger 
than the same line in the spectrum of HD 154090. The same is 
true of the Cm line at 1247 A. However, both stars appear to 
be nitrogen weak (see e.g., the N v wind resonance line around 
1240 A). Therefore on the basis of this evidence, it appears that 
the K Cas sh ould be defined as a ni trogen weak star, rather than 
carbon rich. ICrowther et al.l d2006h found similar results for k 
Cas, citing it as having the 'least nitrogen enriched abundance' 
in their sample as well as the lowest values for the N/C and N/O 
ratios. 



The CNO abundances derived here for the sample of 20 B 
supergiants are compared in Table 7} to values obtained by 
other autho rs ("Trundle et all 120041: iTrundle & Lenn onl l2005[ 
lEvans et alJ l2004a; Crowther et alT 120061; IVenrinT9 95. 199% 
for OBA supergiants. The results from Trundle et al.l (120041) : 



Trundle & Lennon (2005) SMC B supergiants have been com- 
bined to obtain mean CNO abundances based on a sample of 
1 8 stars (but only 1 3 were used for the mean oxygen abundance 
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* = Searle Galactic BSGs 

« = Trundle SMC BSGs 

□ = Crowther Galactic BSGs 

+ = licpolust Galactic sta 




4.Z 

log Tcff 



Fig. 8. Position of the sample of Galactic B supergiants on 
the Hertzsprun g-Russell diagram, a long with other Galactic 

B supergiants ICrowther et alJ (12006). SMC B supergiants 

iTVundle et al. (2004); Trundle & LennonI (l2005h and Galactic O 
stars Repolustet al. (2004). Evolutionary tracks are taken from 
[Maeder & Meynet (2001) and imply 15 Mq < Mevoi < 40 for 
the sample of 20 Galactic B supergiants presented in this work. 
Galact ic B supergiants that are common to both our sample and 
that of ICrowther et al.l (2006) are joined by a dotted red line to 
illustrate the effect of different stellar parameters on a star's pre- 
cise location on the HR diagram. 



since oxygen abundances were not derived for some B2.5 - 5 
stars due to weak, unmeasurable On lines). The data from 
lEvans et al] (l2004ah were purely based on CNO abundances 
derived from OB supergiants so that the results for nebular 
and Hii regions included by the authors for comparison were 
omitted. It is clear from Table |7] that more CNO enrichment 
occurs in stars belonging to the Magellanic Clouds t han Galactic 
stars. This is in accordance with lEvans et al.l(l2004ah . who found 
that OB supergiants in the LMC display a nitrogen enrichment 
that is greater than the nitrogen enrichments in Galactic B 
supergiants. Evans et al. (2004a) conclude that their sample of 
Magellanic Cloud stars show significant nitrogen enrichment 
due to efficient rotational mixing. The CNO abundances show 
no clear trend with effective temperature or sin /. 



4.6. Mass loss rates for B supergiants 

The mass loss rates obtained for this sample of 20 Galactic 
B supergiants are based on matches to the Ha profile and the 
resulting values are listed in Table [8] All of these B supergiants 



typica 



Tor bars arc 



ao 30 40 

Evolutioilai'y Mass 



Fig. 9. Comparison of evolutionary and spectroscopically- 
derived stellar masses for the sample of B supergiants. The dot- 
ted line indicates 1:1 correspondance. 



have mass loss rates ranging between -7.22 < logM < -5.30, 
except the B5 Il/Ib star HD 164353 for which M = 6 x lO *^ 
was derived. The errors quoted in Table [8] reflect the ambiguity 
involved in fitting Ho' 'by eye' (and therefore represent the 
maximum and minimum values of M that fit Ha reasonably) and 
are no greater than a factor of 2. In some cases an upper or lower 
error limit only is quoted where the model fit over- or under- 
estimates the observed Ha profile, meaning that a larger/smaller 
mass loss rate would not be appropriate. CMFGEN fits to the 
Ha profiles of k Cas, HD 190603, HD 14818, HD 190066, 
HD 193183 and HD 164353 are given in Fig. [12] In general, 
good fits are obtained for each star, but several difficulties 
have been encountered in trying to reproduce the observed Ha 
profiles. It is clear that all the observed profiles are asymmetric 
and it is likely that this is caused by the influence of resonant 
line scattering that is too weak to produce a 'P Cygni'-type 
profile so merely results in a slightly asymmetric profile. In 
some stars e.g., HD 213087, it appears to be a redward emission 
component that partly fills in the profile, to such an extent that in 
some stars this red component is visible as a separate emission 
component (e.g., HD 206165) and the Ha profile begins to 
resemble a P Cygni profile (e.g., HD 14818). In the majority 
of stars, the peak/trough of the Ha profile has shifted from the 
line centre as observed in k Cas. This effect is particularly clear 
on comparing the Ha profile of k Cas with that of HD 190603, 
whose peak is much more central resulting in only a slight 
asymmetry to the overall profile. It is also of interest to note that 
CMFGEN predicts a 'bump' in the blueward wing of the Ha 
profile of HD 190603 that is not present in the observed profile; 
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Fig. 10. Comparison of jf^ with luminosity (B5 Ib/II stars 

omitted). Values obtained for the 20 Galactic B supergiants (as- 
terisks) are plotted with those de ri ved for SMC B sup e rgiant s 
(diamonds) iTrundle et aT] (|2004|) : iTrundle & LennonI (|2005|) . 
The dotted line indicates 1 : 1 correspondance. 

Table 5. Calibrations of fundamental parameters by spectral 
type for Galactic B supergiants, based on this work and that of 
[Crowtheret al, (2006) 



Sp. type 




log (L/Lo) 


R,{Re) 




M.(Mo) 


BO la 


28.1 


5.60 


26.9 


2.99 


25 


BO lb 


29.7 


5.66 


23.8 


3.24 


37 


B0.2 la 


26.7 


5.62 


30.4 


3.04 


36 


B0,2 lb 


28.5 


5.65 


27.8 


3.23 


49 


B0,5 la 


24.7 


5.58 


33.8 


2.90 


33 


B0.5 lb 


25.4 


5.58 


32.2 


3.09 


47 


B0,7 la 


23.6 


5.53 


35.1 


2.72 


23 


B0.7 lb 


24.4 


5.51 


33.9 


2.93 


37 


Bl la 


22.0 


5.44 


36.5 


2.41 


12 


Bl lb 


21.7 


5.38 


34.9 


2.67 


22 


B1.5 la 


19.9 


5.44 


44.5 


2.41 


18 


B1.5 lb 


19.3 


5.29 


39.7 


2.50 


19 


B2Ia 


18.3 


5.41 


51.0 


2.32 


19 


B2Ib 


18.1 


5.27 


44.4 


2.46 


21 


B2.5 la 


17.2 


5.39 


56.5 


2.24 


19 


B2.5 lb 


17.6 


5.25 


46.2 


2.43 


22 


B3 la 


16.4 


5.37 


60.4 


2.16 


19 


B3Ib 


17.5 


5.23 


45.5 


2.39 


19 


B4 la 


15.8 


5.34 


63.5 


2.06 


16 


B4Ib 


17.4 


5.18 


43.2 


2.27 


13 


B5 la 


15.7 


5.33 


63.0 


2.03 


15 


B5Ib 


15.2 


5.09 


51.7 


2.11 


13 
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Fig. 11. Comparison of the relative strengths of the N and C lines 
in/f Cas and HD 154090 




Fig. 12. Examples of CMFGEN fits to Ha profiles of k Cas, HD 
190603, HD 14818, HD 190066, HD 193183 and HD 164353, 
along with the mass loss rate adopted for each star. The red, dot- 
ted line represents the CMFGEN model fit to each Ha profile 
and the solid, black line indicates the observed Ha profile. 



a similar phenomenon is observed for HD 193183. 

A small, preliminary investigation into the effects of including 
clumpin g on the morphology of the model Ha profile was un- 
dertaken. IHnHiE&Mnii CSll) assume the winds are clumped 
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Table 8. Stellar wind parameters (M, /3, v^o, Vturb) derived for 
a sample of 20 Galactic B supergiants. The errors given on M 
reflect the errors in fitting each individual Ha profile. 
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Fig. 13. Comparison of CMFGEN derived mass loss rates with 
theoretical mass loss rates predicted by the Vinketal, (200^ 
mass loss prescription. The dotted Une indicates 1:1 correspon- 
dance. 



HDno, Sptype M(10"°Mo yr"') B Voo(km/s) vwb(km/s) 



37128 


BO la 


1 90+"'' 


1.1 


1600 


15 


192660 


BO lb 


5.oo!°:« 


1.3 


1850 


20 


204172 
38771 


B0.2 la 
B0.5 la 


"•-"-0.34 


1.0 
1.1 


1685 
1390 


15 
15 


185859 


B0.5 la 


0.50^1 


1.5 


1830 


20 


213087 


B0.5 lb 




1.5 


1520 


20 


64760 


B0.5 lb 


1 10+'" 


1.0 


1600 


15 


2905 


BC0.7 la 


2 50+" " 


1.5 


850 


20 


13854 


Bl lab 


1 50+"-^ 


1.2 


955 


10 


190066 


Bl lab 




1.5 


1275 


15 


190603 


B1.5 la+ 


2-40!!!:^ 


1.2 


390 


15 


193183 
14818 


B1.5 lb 
B21a 


23+"-" 

-0.00 

1 00+"-^ 


1.5 
1.5 


565 
625 


20 
15 


206165 


B21b 




1.5 


640 


15 


198478 


B2.5 la 




1.2 


550 


20 


42087 


B2.5 lb 


0.50t°° 


1.2 


650 


15 


53138 


B3 la 


0.45^°^ 


1.2 


500 


20 


58350 


B5 la 


70+"-^ 

'"-0.0 


1.0 


320 


20 


164353 


B5 Ib/Il 


0.06^^3 


1.5 


450 


25 


191243 


B5 lb/11 


0.83!":" 


1.0 


550 


20 



with a volume filling factor / and that no inter-clump medium is 
present. The volume filling factor, is defined as: 

f^U + (l.O-U)e-^ (3) 

where /oo is the filling factor at and v^; is the velocity at which 
clumping is 'switched on' in the wind. However, to carry out 
a fair comparison between clumped and homogeneous models 
involves a larger parameter space than merely varying and 
I'd- For example, it is important to check for consistency in the 
atmospheric structure of both models i,e,, that they sample the 
same optical depth in the photosphere and have the same density 
structures, which can include some fine tuning of the velocity 
law and the point at which the TLUSTY hydrostatic structure 
joins the CMFGEN density structure. It is also important to 
check that the same computation options are selected for both 
models to ensure that the density structure is computed using 
the same methods (e,g, the same number of A iterations are 
specified). This is important as computational parameters may 
have been changed for individual models in order to ease 
convergence. For these reasons, a rigorous comparison of 
homogeneous and clumped models will be postponed to a later 
date. 

These mass loss rates have been compared to those p r edicte d 
by the theoretical mass loss prescription of I Vink et al,l (12000 '). 
as shown in Fig, [T3] where the values of T^ff, logL« and M» 
derived in the previous section have been input into the relevant 
mass loss recipes (equations 12 and 13) quoted in the paper 



We find that the M's derived here are in good agreement with 
the I Vink et al,l (l2000h predictions, with discrepancies of a factor 
of 2-3 on average and the maximum discrepancy a factor of 6, 
Trundle et al, (2004); Trundle & Lennon (2005) found that the 
values of Mi,,„t were a factor of five lower than observed mass 
loss rates for early B supergiants, whereas for mid B supergiants 
Mi.ink was a factor of seven higher than observed values. No 
consistent discrepancy is found in our results but generally 
Mvink ^ Miia for BO - Bl supergiants and the reverse is true 
for B2 - B3 supergia nts. Our values of M obtained here were 
compared to those of iCrowther et al,l (12006') . with who we have 
8 stars in common, and the mass loss rates are in very good 
agreement, A co mparison has al so been made to the values 
obtained bv Kudritzki et alj d 19991) . since there are again 8 stars 
common to both data sets (Fig, [T4|, With the exception of e 
Ori, K Cas and H P 206165 , all the B supergiant mass loss rates 
derived by Kudritz ki et al,l (|T999) are smaller than our values by 
typically a factor of up to 5, The values derived for e Ori and k 
Cas are well within the errors of our derived values; however a 
larger discrepancy of a factor of ~ 10 is found for HD 206165, 
Initially this is puzzling since in both cases good fits have been 
obtained to the observed Ha profile of HD 206165 and do not 
suggest such a large discrepancy in M. However, quite different 
stellar parameters have been adopted in terms of T^s iTeff - 
18000 K in our analysis cf, 20000 K from Kudrit zki et al. 
1999), log (L/Lq), and Vk,; more importantlv lKudritzki et al. 
( 1999) adopt a much higher /3 value of 2,5 compared to 1.5 in 
this work. 
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-6,4 -6.2 -6,0 -5.8 

log mass loss (Scarlc) 



Fig. 14. Comparison of CMFGEN derived M with those of 
[Kudritzki et al.. (.1999.) for 8 stars common to both samples 



4.7. The Wind-Luminosity-Momentum Relation 

The concept of a Wind-Luminosit y-Momentum Relation (here- 
after WLR) was first proposed bv lKudritzki et al.l (Il995h . using 
the prediction from the theory of radiatively driven winds that 
there is a strong dependence of the total mechanical momen- 
tum flow Mvoa of the stellar wind on stellar luminosity (e.g., 
ICastor et al.lll975h . which can be described as 



Mvoo oc RjL°'ff 



(4) 



The importance of the WLR lies in its potential as an extra- 
galactic distance indicator provided that it is reliably calibrated. 
The proportionali ty shown i n Equ ationlHwas first confirmed ob- 
servationally by Puis et al.l ( 11996! ?) for a sample of Galactic and 
Magellanic Cloud O stars with 5.5 < log L» < 6.5. For log L» < 
5.5, a linear fit was not possible, demonstrating the d epen- 
dence of the WLR on spectral type. iKudritzki et al.l (1 19991) then 
showed that a linear fit to the WLR was also pos sible for galac- 
tic BA supergiants. Since then many authors ( Repolustetal 
2004; Markova et al. 2004; Massey et al. 2004; Masse v et al 



20051; [Trundle et al.il2004l; iTrundle & Lennonl l2005) have pub- 



lished values for wind momenta when deriving fundamental pa- 
rameters for sets of OBA stars using 



D„ 



= Mv^R 



,0.5 



(5) 



where is in solar radii, M in g/s and Voo in cm/s to give D„,om 
in units of cgs. Assuming a WLR of the form 

logD„,o„, = logA) + x\og (L/Lq) (6) 



a linear regression can be used to constrain the coefficients x 
and log Dq. The reciprocal of x can be thought of as the effective 
exponent a^ff (see Equation |4]i. Applying a linear regression 
to our data gives log Do - 19.76 and x = 1.61. Looking at 
Fig. [TSj it can be seen that the hotter spectral types i.e. O stars 
have a steeper WLR than cooler B spectral types. This is to 
be expected if, as proposed by IVink et al] (Il999l) . Feii and 
Fe III lines are responsible for driving the subsonic part of the 
wind, corresponding to lower ionisation stages. Fig. [15] also 
illustrates the effect of metallicity on the WLR, with the more 
metal poor environment of the Magellanic Cloud hosting stars 
with lower values of D,,,,,,,,. This effect is particularly notice- 
able between Galactic and SMC B supergiants. Our values 
of Dmom are compar ed to those pred icted by the theoretical 
WLR prescription of ' Vink et all (1200 0), using the parameters 
derived for our sample of stars in this chapter. It is found that 
the observational values are greater than predicted values for 
BO - B0.7 supergiants (except for the B0.2 la star HD 204172) 
where T^ff > 23 000 K, which is expected for the hotter side 
of the bi-stability jump. Conversely, B 1 - B5 supergiants have 
smaller, observed values of Dmom compared to predicted values. 
This is caused by a combi nation of several effects. Firstly, just 
as the predicted lVink et akl ([2000.) mass loss rates are a factor of 
5 larger on the cooler side of the bi-stability jump, (reff ~ 23 000 
K), the predicted wind momenta will be greater for Bl - B5 
supergiants, causing a larger discrepancy between observed and 
theoretical values of Z)mom- In addition to this, many Bl - B5 
supergiants have Ho- profiles in absorption, making it harder 
to constrain a 'true' observed m a ss los s rate. Similar results 
are found by Trundle & Lennon (2005) for their sample of 
SMC B supergiants. Repolustetal. (2004) suggested that the 
inclusion of clumping in the derivation of mass loss rates may 
help to alleviate the existing discrepancies between observed 
and theoretical wind momenta, which also exist for O stars as 
well as B stars. Clumping would reduce the mass loss rate and 
consequently lower wind momenta, although the precise amount 
of clum ping present in O B star winds remains uncertain. Very 
recently TPuls et al.) (120061) attempted to derive better constraints 
on the clumping factor in hot star winds through a combined 
optical, infra-red and radio analysis of the wind. They found 
that use of clumped mass loss rates did produce much better 
agreement between observed and theoretical wind momenta, 
since for O stars the observed wind momenta were originally 
higher than the theoretical ones, but the inclusion of clumping 
reduced the value of M and consequently D„om- Unfortunately, 
for the case of Galactic B supergiants this may only work for 
BO - BO. 7 supergiants; the use of a lower, clumped M would 
not resolve the discrepancy for Bl - B5 supergiants. This may 
indicate a fundamental difference between the structure and 
inhomogeneity of O and B star winds. At present, it is not 
possible to obtain a reliable calibration of the WLR, until the 
problems associated with its dependence on luminosity and 
metallicity, as well as the effect of clumping on observed mass 
loss rates, are resolved. 



5. Testing the UV predictions of CIVIFGEN 

The next step of our investigation was to examine if the 
CMFGEN models presented in the previous section would 
also provide a good fit to the UV silicon lines, thus confirming 
that Tetf diagnostics at both optical and UV wavelengths 
implied the same value of T^ff for each star. An example of a 
CMFGEN comparison to Teff-sensitive silicon lines Sin /11265, 
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Fig. 15. The Wind-Luminosity Momentum Relation for OBA 
stars. Note the dependence of the WLR on spectral type and 
metallicity. The theoretical WLR predicted by Vink et al.l (l2000l) 
is calculated for our sample of Galactic B supergiants and is rep- 
resented by the orange line. 



Sim /1 1294 and ^1299 (as noted by iMassal [19891) is giyen in 
Fig. [16] for the B2 la star HD 14818 and the corresponding 
final CMFGEN model with Teff = 18 000 K, L = 2.5 xlO^ 
L/Lq and M - 1.1 x 10"^ Mq yr Two other models with 
Teff = 17500 K, L = 2.4 xlO^ L/Lq, M = 1.2 x 10"^ Mq yr"' 
and Teff = 18 500 K, L = 2.5 xlO^ L/Lq, M = 1.8 x 10"^ 
Mq yr"' respectiyely are also shown to demonstrate the effects 
of changing Teff on these lines (the slight differences in the 
luminosity and mass loss of these models will not significantly 
affect the silicon lines). A direct comparison of the obseryed and 
model Sin /il265 line profiles is difficult since the continuum 
is raised about this line, but it is apparent that the model 
produces an asymmetric profile (whereas the obseryed profile 
is symmetric) shifted by about 1 A blue-ward relatiye to the 
obseryed line profile centre. The model profile is also much 
broader than the obseryed profile and yarying T^ff by + 500 K 
has no significant eff'ect on this line. In the case of Sin /il309, 
the model line profile is more narrow and shallow than the 
obseryed one. Changes in Teff are more apparent on this line, 
though still make no significant improyement to the oyerall line 
fit. For Sim A1294 and /1 1299, the model line profiles are again 
asymmetric, unlike the obseryed profiles, and the blue wings of 
these lines are oyerestimated whilst the absorption troughs are 
underestimated. In fact, both obseryed Si m lines appear to show 
some eyidence of broadening due to the stellar wind despite 
being photospheric, which is also eyident in the model profiles 
in the form of asymmetry. Additionally change in Tetf appears 
to haye no aff'ect on these lines; howeyer the higher yalue of 




Si III 1294 Si III 1299 




Wavelength 



Fig. 16. CMFGEN model fit to the lUE spectrum of HD 14818 
(B2 la), focusing on the UV silicon T^ff diagnostics: Sin /11265 
and AUQ9, Sim ^1294, ^1299 and ^1417. The 'best fit' model 
deriyed from the optical has Teff = 18 000 K, L = 2.5 xlO^ L/Lq 
and M - 1.1 x 10"* Mq yr"' (dashed red line). Others models 
hayeTeff = 17 500 K,L = 2.4x10^ L/Lq,M= 1.2x10"'' Mq yr"' 
(dotted green line) and T^ff = 18 500 K, L = 2.5 xlO^ L/Lq, M 



1.8 X lO"*" Mo yr"' (dot-dashed blue line). 



M for the T^s = 18 500 K model (blue line) produces a deeper 
absorption trough for the profile. To conclude, yarying T^ff and 
eyen M has a small affect on these lines, but will not succeed 
in reproducing the obseryed lines accurately, with the correct 
broadness and symmetry. 



5.1. The optical/UV discrepancy in wind lines 

A large number of UV lines are also strongly affected by mass 
loss from the wind, so it is also of interest to inyestigate whether 
the yalues of M deriyed from Ha in 94.6l succeed in reproducing 
the UV wind resonance lines accurately. This is not the first 
time that modelling of hot stars has been extended to the UV 
and matching the P Cygni profiles obseryed there. In the last 
couple of years, seyeral authors haye begun to consider both 
optical and UV stellar properties when denying fundamental 
parameters (e.g., Eyans et al. 2004a; Crowther et al. 2006) and 
iBouret et al.l (12005 ) analysed lUE and FUSE spectra of two 
Galactic 04 stars with CMFGEN and TLUSTY, presenting one 
of the first analyses based exclusiyely on UV diagnostics that 
also uses these particular stellar atmosphere codes. It is eyident 
that an optical analysis proyides a much easier way of obtaining 
stellar parameters, where diagnostics for e.g. T^ff and luminosity 
are readily ayailable and only depend on abundance, T^ff and/or 
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luminosity. On the other hand, the task of identifying suitable 
diagnostic lines is less straightforward, since many UV lines 
will be sensitive to mass loss as well as T^ff, abundance and in 
some cases Vturb- 

An example of a CMFGEN fit to the lUE spectrum of 
HD 190603 (B1.5 Ia+) is shown in Fig. [TT] Since the mass 
loss rate has already been constrained from fits to the Ha 
profile, it is interesting to see whether the derived value of M 
is confirmed by reasonable fits to the UV P Cygni profiles, 
provided that a reasonable model fit to the Ha profile has 
already been achieved. Looking at the case of HD 190603 
shown in Fig. [ITl the fit to the observed Ha profile is good. 
However, it is clear that CMFGEN does not reproduce any 
of the observed P Cygni profiles accurately, implying that a 
different value of M would be appropriate for the UV. The 
model fails to produce sufficient high velocity absorption in 
the UV wind resonance lines, to the extent that the predicted 
C IV AA 1548.2, 1550.8 line is only present as a photospheric line 
with no evidence of wind contamination. Nv is not seen as a 
P Cygni profile in this star, but the model does not even produce 
a distinct, weak photospheric line at 1238 A. However, better 
fits are achieved at lower ionisation: Cii AA 1335.66, 1335.71; 
Si IV AA 1393.8, 1402.8 and Aim AA 1854.7, 1862.8. The 
observed Cii AA 1335.66, 1335.71 line is saturated but the 
model produces an unsaturated line, which suggests that either 
a model with a higher value of M is required or the model 
ionisation is incorrect. Adopting a higher value for M though 
would worsen the effect of the model overestimating the red 
wings of the Siiv AA 1393.8, 1402.8 doublet. It would have a 
more positive effect on the Aim AA 1854.7, 1862.8 line, since 
the observed blueward doublet is beginning to saturate but the 
model blueward doublet is clearly unsaturated, again supporting 
a higher mass loss rate. The CMFGEN fit to Ha would worsen 
if a higher value of M was adopted, illustrating the discrepancy 
between the mass loss rates implied from the optical and UV. 
It is also noticeable when comparing the observed and model 
Si IV P Cygni profiles that the model doublet components are 
narrower than observed. As in the case of Civ, this is due 
to the model predicting to little absorption at high veloci- 
ties. For HD 190603, these problems arise in spite of the fact 
that the value adopted for M provides a good fit to the Ha profile. 

An example of a better CMFGEN fit to the UV wind resonance 
lines is given in Fig. [18] for the B2 la star HD 14818. The 
observed Ha profile displays a P Cygni profile, which has not 
been successfully reproduced by the model (as discussed in 
34.61) . Despite this, very good fits have been obtained to Si iv 
and Aim in comparison to those obtained for HD 190603, 
though again a lack of high velocity absorption causes the 
model to under-estimate the broadness of the absorption trough 
for Si IV. However, the same failure occurs in reproducing the 
P Cygni profile of Civ line, whilst Nv shows no evidence of 
wind contamination. The fit to C ii is reasonable, although the 
model predicts too much redward emission and as a result does 
not match the redward side of the absorption trough. Conversely 
an example of a worse fit than either of the previous cases 
is shown in Fig. [19] for HD 53138 (B3 la). Its observed Ha 
profile is in absorption but shows a small amount of red-ward 
emission and is reasonably well matched by CMFGEN. On 
the other hand, the UV P Cygni profiles are in general poorly 
matched by the model, with none of the five wind line profiles 
being well reproduced. The same problems seen for HD 190603 
and HD 14818 in matching N v, C iv and Si iv also occur here. 




C IV 1640-2, 1&50.S 




H alpha 6555 

J 







Fig. 17. CMFGEN fit to the lUE spectrum to the N v, C iv. Si iv. 
Aim and Cii wind resonance lines of HD 190603 (B1.5 la-i-). 
Note that a good fit to Ha does not guarantee the same mass loss 
rate will provide a good fit to the UV P Cygni profiles. Model 
pai-ameters ai-e T^s = 19500 K, log (L/Lq) = 5.41 and M = 2.4 
xlO"*" Mo. 




Fig. 18. CMFGEN fit to the lUE spectrum to the N v, C iv. Si iv. 
Aim and Cii wind resonance lines of HD 14818 (B2 la). Even 
though the fit to Ha is not perfect, a reasonable fit is made to the 
UV P Cygni profiles, particularly Si iv and Al iii. Model param- 
eters are T^ff = 18 000 K, log (L/Lq) = 5.40 and M = 1.0 xlO"^ 
Me. 



The red-ward emission in Cii is grossly over-estimated and 
the model produces an asymmetric Aim profile that is not 
observed. In both cases the model predicts saturated lines when 
the observed profiles are not saturated (though C ii is beginning 
to saturate a little). The high velocity absorption in Aim is 
over-estimated to the extent that it predicts saturation to occur 
at a higher velocity than observed. It is therefore clear from Fig. 
[TT] Fig. [18] and Fig. [19] that a discrepancy exists between the 
value of M required to fit the Ha and UV wind resonance lines 
(hereafter referred to as the opticalfUV discrepancy). 



S. C. Searle et al.: Fundamental parameters of Galactic B supergiants 



19 




Fig. 19. CMFGEN fit to the lUE spectrum to the Nv, Civ, 
Si IV, Aim and Cii wind resonance lines of HD 53138 (B3 la). 
Although a good fit has been made to Ha with the adopted 
mass loss rate, CMFGEN does not reproduce the observed UV 
P Cygni profiles well. Model parameters are Tejf - 16500 K, 
log (L/Lo) = 5.30 and M = 4.5 xlO"^ Mq. 



In general, CMFGEN only succeeds in matching the C iv line 
when it is saturated in early B supergiants, at which point it is 
no longer sensitive to Teff and M so a reliable fit cannot be ob- 
tained as altering these parameters will have no affect on the 
model line profile. Otherwise, CMFGEN manages to reproduce 
most of the observed P Cygni profile for C n, Al iii and Si iv, 
but fails to produce enough high velocity absorption to repro- 
duce the full extent of the observed absorption trough. As a re- 
sult, the model often under-estimates the blueward absorption 
as well as over-estimating the redward emission, especially in 
the case of Si iv. This can sometimes lead to the model giving 
an asymmetry to the P Cygni profile that is certainly not ob- 
served in the spectrum. Additionally, CMFGEN never succeeds 
in producing the Nv P Cygni profile when present in BO - Bl 
supergiants and even when a weak, photospheric profile is ob- 
served, the model fails to produce a discernible spectral line 
at the correct wavelength for N v. In the hotter B supergiants, 
the model grossly underestimates the photospheric Al iii and C ii 
lines. However when the same resonance lines are seen as P 
Cygni profiles, the model has a tendency to reproduce them as 
saturated when they are observed to be unsaturated. All these 
discrepancies suggest that the problem lies within the predicted 
ionisation structure of the models. CMFGEN fits to the overall 
lUE spectra of 10 BO - B5 supergiants are available as online 
material (Fig.s|9l [T0l[TTT i. Very similar problems in matching the 
UV P Cygni profiles have als o be en encountered by Evans et al. 
( l2004al) and iCrowther et all (|2006|) when modelling O and early 
B supergiants with CMFGEN. 

5.2. Modelling the UV exclusively 

The CMFGEN models examined in the last section demonstrate 
a clear discrepancy between Mua and the value of M implied 
by the P Cygni profiles of the wind resonance lines. It is 
hardly surprising that they are unsuccessful in reproducing the 
observed UV wind diagnostics accurately. In this section, the 



possibility of modelling a star solely from its UV spectra will 
be investigated (ignoring any prior knowledge of values of 
parameters from the optical) to see if the UV can be reproduced 
more accurately. In order to do this, we must first identify 
suitable UV diagnostic lines by which values of Teff, log (L/Lq), 
M, Voo, /3 and abundances could be constrained. 

Looking back to the problems mentioned in the previous sec- 
tion, one potential difficulty is immediately apparent. CMFGEN 
is unable to reproduce the Civ line accurately, which makes 
it hard to constrain Vm and /3 from this line. Suitable UV Teff 
diagnostics also need to be found besides the photospheric Si ii 
and Sim lines discussed in ^ Siiv AA 1393.8, 1402.8 could 
be a good candidate but it is also very sensitive to luminosity 
and mass loss; moreover it is often saturated, reducing its 
sensitivity to both parameters, and CMFGEN rarely reproduces 
it accurately. Other potential T^g diagnostics are Al iii and C ii 
which also show some sensitivity to mass loss and are therefore 
not ideal. Another possible Tgff diagnostic is the photospheric 
Si II 1526.7, 1533.4 A line, but CMFGEN does not model these 
lines well either, often completely failing to reproduce the 
blue-ward part of the doublet. More importantly, the 1533.4 A 
doublet becomes blended with C iv AA 1548.2, 1550.8 A at high 
value of Voc- At this stage, we have no photospheric lines to 
use as reliable Teft diagnostics, since they are not well matched 
by CMFGEN. The best we can do is look at the UV lines best 
reproduced by CMFGEN (i.e., Siiv, Aim and Cii) and analyse 
their sensitivity to the main stellar parameters. 

In practice, another major problem materialises. It is difficult to 
disentangle the effects of Tgff and M on Si iv, Al iii and C ii, plus 
they are often too saturated to be sensitive enough to these pa- 
rameters. When Si iv is not observed to be saturated, CMFGEN 
still predicts a saturated profile that is virtually insensitive to Teff 
and M, making it difficult to use as a T^^fi and M diagnostic. In 
fact, the lack of a significant difference between model P Cygni 
profiles when varying mass loss presents a serious obstacle 
to any attempt to derive parameters from the UV, as we will 
now show. For BO-Bl supergiants, the model often produces 
a saturated C iv P Cygni profile and over-estimates the Si iv P 
Cygni profile. It may appear logical that adopting a model with 
a lower mass loss rate would provide a better fit to the observed 
Civ and Siiv lines. However, the lack of sensitivity of this line 
to mass loss becomes apparent when the M adopted by the 
model is altered. This is illustrated in Fig. |20l where it can be 
seen that lowering the value of M from 5.0 x 10"^ to 2.6 x 10"* 
has no affect on the wind resonance lines (implying that they 
are still optically thick), despite producing model Ha profiles 
in emission and absorption respectively (note that the broad 
feature seen in the model between 1242 - 1247 A is not N v but 
C III, which interestingly enough does show some sensitivity to 
mass loss). It could still be argued that a larger decrease in mass 
loss is required to fit these lines. However Fig. |2T| disproves 
this idea as yet again no difference is seen between P Cygni 
profiles for models with M = 6 x 10"^ M© yr"' (red dashed 
line) and = 1.8 x 10"^ M© yr"' (blue dotted line) respectively. 
This is in spite of the fact that this difference in mass loss again 
results in model Ho' profiles in emission and absorption, as 
well as having a significant difference on the amplitude of Ly 
a (1216 A). HD 164353 presents an interesting case study for 
how CMFGEN deals with the ionisation in the stellar wind, 
as it is a B5 Ib/II star that possesses a very weak wind with 
M = 6 X 10"^ Mq yr"' and can be thought of as a star with 
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Fig. 20. Comparison of UV wind resonance lines of HD 192660 pj ^i. Comparison of UV wind resonance lines of HD 164353 
for models with M = 5 x lO"'' (red dashe^^^^^^ and 2.6 x ^^^^j^ ^.^j^ M = 6 x lO ** M^, yr"' (red dashed line) and 



10" Ma (blue dotted line) for HD 192660 (BO lb). 



= 1.8 X 10 Mq yr (blue dotted Une) respectively 



negligible mass loss and stellar wind contamination. This is 
confirmed by looking at the observed C ii and Al iii resonance 
lines (Fig. \2l[ . which show some asymmetric broadening. 
However the model predicts strongly saturated profiles for both 
lines despite the low mass loss rate adopted for the model, 
again suggesting that the predicted ionisation structure is at 
fault. This highlights another significant problem that, given 
their lack of sensitivity to significant changes in mass loss, the 
Civ and Siiv P Cygni profiles would not make suitable mass 
loss diagnostics. It appears that the root of the problem lies in 
CMFGEN predicting ionisation fractions for C ii and Al m that 
are too high, resulting in a large optical depth that produces too 
many absorbers at too high a velocity. The observed profiles 
on the other hand show us that absorption is only occuiTing 
around the rest velocity of the line. The model over-estimation 
of C II and Al iii may therefore only be resolvable by lowering 
the ionisation fraction of these two elements and cannot be 
resolved by altering the mass loss rate of the model in question. 
From this, we conclude that the UV wind resonance lines are 
not suitable candidates for deriving T^^fi and M. Even if the 
ionisation structure was correctly predicted, more diagnostic 
lines would be required to determine all the necessary stellar 
parameters other than T^ff and M, as well as ensuring that an 
accurate analysis had been carried out. 

In addition to the wind resonance lines, the UV subordinate 
lines can potentially be used to provide additional constraints 
on the mass loss adopted for the model. An example is the 
Si IV 1122, 1128 A line in the FUV, whose upper energy level 
is coincident with the lower energy level of Siiv 1400 A. 



This means that if the model over-populates the lower level of 
Siiv 1400 A, the upper level of Siiv 1122, 1128 A will also 
be over-populated, pushing the line into emission when it is 
observed to be in absorption. If this predicted line is seen to be 
in emission in a model when the observed line is in absorption, 
this is a direct indication that the adopted mass loss rate of the 
model is too high. There are no examples of Si iv 1 128 A being 
in emission in the models used for the sample of 20 Galactic B 
supergiants, so this would suggest that the adopted mass loss 
rates are within reason. This also means that we could not have 
used this line as a mass loss diagnostic in this analysis. 

It is possible to provide alternative perspectives on the UV be- 
haviour by comparing the ionisation stages present in any given 
star, rather than trying to reproduce individual line profiles. Fig. 
|22] shows the predicted ionisation structure against w at four 
different Teft ; 27 500 K, 23 500 K, 18 000 K and 15 000 K for the 
six ions (N v, C iv. Si iv. Si iii, Al iii and C ii). CMFGEN predicts 
that Si IV will be dominant as expected for Tgif = 27 500 K, but 
shows very low levels of N v and C iv. This is hardly surprising 
since it explains the complete absence of a N v P Cygni profile 
(when present observationally), as well as the difficulties in 
generating a P Cygni profile for C iv when it is unsaturated. It 
is also interesting to note that the levels of ionisation drop off 
rapidly in the mod el as w increa s es, co ntradicting the empirical 
determinations of iPrinia et akl (l2005b where winds became 
more highly ionised at high w. This is direct evidence of the 
model failing to generate enough high-velocity absorption to 
sustain the same level of ionisation further out in the wind. 
This is the reason for the 'narrowness' of the model Civ and 
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Fig. 22. CMFGEN predicted ionisation structure at different T^fi, 
plotted against normalised veocity w. The ions are colour-coded 
as follows: N v = green, C iv = red, Si iv = dark blue. Si iii = light 
blue, Al m = purple and C ii = yellow. Si iv is predicted to be 
dominant for B0-B2 supergiants (30000/: < T^ff < ISOOO/T), 
after which Si iii, Al iii and C n take over as the dominant ions in 
the wind for T^ff < 18 000 K. 



Si IV P Cygni profiles compared to the broad absorption troughs 
of the observed P Cygni profiles. If the model cannot sustain 
enough ionisation in the inner and outer parts of the wind, then 
it will be unable to fully reproduce the blue-ward part of the 
profile. CMFGEN predicts Siiv to be dominant down to Teff 
= 18 000 K, at which point Sim and Aim take over as the 
dominant ions in the wind. At this Tgff, C ii has also increased in 
strength, becoming a dominant ion at Teff = 15 000 K. Whereas 
this approach has provided us with valuable insight into why 
CMFGEN struggles to predict the P Cygni profiles correctly, it 
too fails to provide us with an alternative means of constraining 
parameters due to the incorrectly-predicted ionisation structure. 

Given all these problems with CMFGEN mismatching the 
observed UV P Cygni profiles, a investigation into the effects 
of clumping on these lines would not be worthwhile at present. 
In addition, mass loss in the UV is only sensitive to p, rather 
than as in Ha and radio-dominated regions of the wind, so 
it is not a particularly sensitive indicator of clumping. First the 
models need to predict the correction ionisation structure for B 
supergiants. Secondly, the problems associated with investigat- 
ing the effects of clumping on Ha (as discussed in § 14.61 1 need 
to be sorted as Ha is an important diagnostic of clumping and 
can provide important insight into its behaviour, which would 
aid a subsequent analysis of clumping in the UV. Furthermore, 
in comparison to O stars which possess strong indicators of 



UV clumpiiig e.g. P v /li 1 1 18, 1128 (see lFuUerton et al.ll2006l: 
iBouret et al.l 120051) . B supergiants do not possess an equally 
convenient UV diagnostic. A tentative examination of the effect 
of clumping on the UV profiles has shown that it does not 
improve the fits to the observed P Cygni profiles, as expected, 
but can alleviate the over-estimation of emission seen in the 
red-ward part of the model Si m and Al m profiles. In the case of 
the photospheric Sim lines around ~ 1300 A, some broadening 
of these lines due to the stellar wind is seen observationally, 
and the models also show some sensitivity to clumping in these 
lines. The models exaggerate the effect of the stellar wind on 
these lines by producing slightly asymmetric profiles, but the 
inclusion of clumping can help to lessen this asymmetry. This 
is logical since the inclusion of clumping will increase the 
wind density locally, providing more absorption at the point at 
which the line forms in the wind, helping to reduce the excess 
red-ward emission seen in many of the model P Cygni profiles. 
The inclusion of clumping will have no affect on saturated lines 
in the model as they are no longer sensitive to density changes 
in the wind. All the afore-mentioned issues associated with 
investigating clumping effects need to be addressed before any 
truly meaningful analysis of clumping in the UV can be caiTied 
out. 



6. Conclusions 

A quantitative study of the optical and UV properties of BO - B5 
la, lab, Ib/II supergiants has been carried out, usi ng the nUTE, 
line-blanketed stellar atmosphere code of iHillier & Miller] 
(1998). A revised B supergiant Teff scale (derived using a stellar 
atmosphere code that includes the effects of line blanketing) has 
been presented, giving a range of 14.5 000 K < Tgff < 30000 K 
for these stars. This scale shows a drop of up to lOOOOK from 
BO la/b to B 1 lab and a difference of up to 2 500 K between la 
and lb stars. It also shows that on average the effect of including 
line blanketing in the model produces a modest reduction of up 
to 1 000 K for BO - BO. 7 and B3 - B5 supergiants, whereas a 
larger reduction of up to 3 000 K is seen for B 1 - B2 supergiants 
(see Table |4]i. The 20 Galactic B supergiants also displayed a 
range of 2.1 < log g < 3.4 in surfac e gravi ty. These results, 
together with those of ICrowther et al.l (l2006l) . have been used 
to construct a new set of averaged fundamental parameters 
for BO - B5 supergiants, according to spectral type. Mass 
loss rates derived from Ha proved B supergiant winds to be 
generally weaker than those of O supergiants (as expected 
since they are lower-luminosity objects) with M ranging from 
-7.22 < logM < -5.30. All 20 B supergiants also shown signs 
of CNO processing, with the largest nitrogen enrichments being 
seen for B1-B2 supergiants. Evidence for a mass discrepancy 
is found between estimates of Mspec and Mevoi, with the largest 
differences peaking at a value of log (L/Lq) ~ 5.4. 

A Wind-Momentum-Luminosity relation has also been derived 
for our sample, which is lower in val u e for B 1 - B5 supergiants 
than that predicted by IVink et all (12000 ^. but greater than 
predicted values for BO - BO. 7 supergiants. For this reason 
it is not possible to reconcile this difference in observed and 
theoretical WLRs over the whole B supergiant spectral range 
by adopting clumped M as is the case for O stars. A severe 
problem exists in the form of the optical-UV discrepancy, 
where the model fails to reproduce some of the P Cygni profiles 
accurately. This highlights a failure in the model to generate 
enough high-velocity absorption to succeed in reproducing the 
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observed P Cygni profile and more crucially highlights that the 
models are not predicting the correct ionisation structure. Given 
that B supergiants, along with other massive stars, have their 
peak flux in the UV, it is imperative that this discrepancy is re- 
solved if we are to have confidence that fundamental parameters 
derived by this method are a true representation of the star's 
properties. Furthermore it underlines the incompleteness of our 
current understanding of the physics of massive star winds and 
the necessity to review the standard model. A more thorough 
analysis of the ionisation structure of early B supergiant winds 
will be presented in Paper II. 
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Appendix A: Error Analysis 

In this section, the errors affecting each derived parameter are 
discussed. The error on Teff is estimated from the quality of the 
CMFGEN model fit to the diagnostic silicon, helium and magne- 
sium lines and therefore represents the range in Tgff over which 
a satisfactory fit to the observed spectrum of the star could be 
obtained. Luminosity is primarily constrained through dered- 
dening the observed spectra with respect to the model spectral 
distribution, its eiTor depends on AMy, whose eiTors are esti- 
mated from dereddening the observed spectrum with respect to 
the model spectral energy distribution. Alog L, also depends on 
AMboi, since Mboi = My + B.C., therefore AlogL» is calculated 
as 

2.5AMboi .. 
AlogL, = logL, (A.l) 

The error on /?, depends on the square root of the sum of 
(ATeff)^ and (AL,)^, with AL» having the greatest influence on 
AR^. For 15 out of the 20 B supergiants, AR^ is within 10% 
of the absolute value of those stars with larger AR» are 
discussed separately in this section. The error on log g for 16 
of the 20 stars star is estimated to be 0.25 dex, based on the 
accuracy of line fits and the effect of AT^ff in determining log g. 
For four of the stars in our sample, we adopted A log g = 0.38 
dex. In the cases of HD 190603 and HD 14818, this was due 
to asymmetric nature of the Hy and H6 profiles, whereas for 
HD 64760 and HD 13854 it reffected the larger eiTor in T^b of 
2000 K (cf. to 500 - 1500 K for other stars in the sample. The 
resulting uncertainty on the spectroscopic mass, Mspec, due to 
errors in constraining log g and range from 0.08 < AM, < 
14.52. In comparison, AMen,i is typically 5 Moas shown in 
Fig. [8] which demonstrates the effect of assuming different val- 
ues of Teff and log {L/Lq) on a star's position in the HR diagram. 

Determining the eiTor in constraining the mass loss rate is 
more complicated, since it depends on both AR^ and the error 
in fitting the Ha profile by varying M and /3. However, the 
errors incurred from uncertainties in deriving are negligible 
compared to those arising from fitting the Ha profile, so we are 
justified in defining AM as solely the eiTor in fitting the Ha 
profile, accounting for the degeneracy in varying /3 to fit Ha 
profiles in emission. Values of Vk, are taken from SEI analysis of 
UV wind resonance lines (the result of which will be presented 
in Paper II) and are accurate to + 50 km/s. The values for Vmrb 
are constrained with an uncertainty of ± 5 km/s, as dictated by 



sensitivity of fitting the Si iii lines by eye. 

Some uncertainties exist in our analysis that waiTant further 
discussion. Although for the majority of stars it was possible to 
constrain T^ff within ± 1 000 K, this was not possible for the 
stars HD 64760 and HD 13854. HD 64760 is a rapid rotator 
and the large width of its spectral lines makes it harder to make 
an accurate distinction between model fits whose T^ff differ by 
1000 K, resulting in A Teff = + 2000 K for this star. In the 
case of HD 13854, if the adopted T^ff of 20000 K is increased 
to 22 000 K (keeping the same luminosity) then a much better 
fit is made to the silicon lines (i.e., Siiv 4089 and Sim 4552, 
4568 and 4575) but at the expense of grossly over-estimating 
the hydrogen and helium lines (i.e., Hy6, Hy, H6, He: 4121, 
4144, 4387 and 4471). Normally in our analysis, the fitting of 
silicon lines would be given priority, but given the weakness of 
Si IV 4089 at this spectral type (Bl lab), it is reasonable to assign 
greater importance to fitting the helium lines. Furthermore, 
adopting a model that only fits the silicon lines well and 
largely over-estimates the hydrogen and helium lines will give 
a misleading indication of the value of T^ff. HD 13854 also 
has quite a large error in My, which consequently propagates 
into significant uncertainties in log (L/Lq) and, combined with 
a larger Areff, leads to a very large AR,. This arises from a 
noisy lUE SWP spectrum and the absence of a LWR spectrum 
leading to a large dereddening error; the same is true for the 
considerable eiTors on the values of My and /?, obtained for 
HD 204172. Low quality lUE spectra generating higher A 
E(B-V) also explain the Alog (L/Lq) found for HD 14818 and 
HD 206165. However, for the B5 II (lb) star HD 164353, it is 
the value of My = -4.2 that poses a problem; CMFGEN simply 
fails to calculate a succesfully-converged model at the required 
luminosity. This explains the large values of Alog (L/Lq), 
A/?, and AMy. The adopted value of My has been independently 
confirmed by several different sources in the literature so we 
believe it to be correct. Furthermore, four stars (HD 37128, 
HD192660, HD198478 and HD 42087) have larger errors in 
the observed value of My than the value of My derived from 
dereddening, so in practise the quoted value of AMy could be 
up to 0.2 mag larger. 
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Spectra of B la supergiants: 4050 - 4250 A 
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Fig. 1. CMFGEN model fits (4050 - 4250 A) to the optical spectra of 10 B la supergiants, with the T^s, luminosity and CNO 
diagnostic lines marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 




Fig. 2. CMFGEN model fits (4250 - 4450 A) to the optical spectra of 10 B la supergiants, with the Teff, luminosity and CNO 
diagnostic Unes marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 
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Spectra of B la supergiants: 4450 - 4650 A 
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Fig. 3. CMFGEN model fits (4450 - 4650 A) to the optical spectra of 10 B la supergiants, with the T^g, luminosity and CNO 
diagnostic lines marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 
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Spectra of B lb supergiants: 4050 - 4S50 A 
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Fig. 4. CMFGEN model fits (4050 - 4250 A) to the optical spectra of 10 B lb supergiants, with the Tes, luminosity and CNO 
diagnostic lines marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 
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Spectra of B lb supergiants: 4350 - 4450 A 
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Fig. 5. CMFGEN model fits (4250 - 4450 A) to the optical spectra of 10 B lb supergiants, with the T^s, luminosity and CNO 
diagnostic lines marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 
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Fig. 6. CMFGEN model fits (4450 - 4650 A) to the optical spectra of 10 B lb supergiants, with the T^s, luminosity and CNO 
diagnostic lines marked as shown. Optical spectrum is in black, CMFGEN model fit is shown in red. 
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Fig. 7. TLUSTY model fits to the Hy profile of all 20 supergiants. Optical spectrum is represented by a solid, black line; TLUSTY 
model fit is shown as a dotted black Une. 
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Fig. 8. TLUSTY model fits to the H6 profile of all 20 supergiants. Optical spectrum is represented by a solid, black line; TLUSTY 
model fit is shown as a dotted black hne. 
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Fig. 9. CMFGEN model fit to the lUE spectra of 10 BO - B5 supergiants (1230 A - 1480 A). The solid red line represents the model 
fits whereas the solid black line is the lUE spectrum. 
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Fig. 10. CMFGEN model fit to the lUE spectra of 10 BO - B5 supergiants (1480 A - 1680 A). The solid red Hne represents the 
model fits whereas the solid black line is the lUE spectrum. 
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Fig. 11. CMFGEN model fit to the lUE spectra of 10 BO - B5 supergiants (1680 A - 1880 A). The solid red Hne represents the 
model fits whereas the solid black line is the lUE spectrum. 



